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c3 1 ABSTRACT: 



In this report we consider the numerical simulations at finite temperature using lattice QCD 
data for the computation of the thermodynamical quantities including the pressure, energy 
density and the entropy density. These physical quantities can be related to the equation 
of state for quarks and gluons. We shall apply the lattice data to the evaluation of the 
specific structure of the gluon and quark condensates at finite temperature in relation to the 
deconfinement and chiral phase transitions. Finally we mention the quantum nature of the 
phases at lower temperatures. 
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I. Introduction: Thermodynamics of Strong Interactions 



We start our considerations after this opening discussion by describing the numerical 
calculations on the lattice using the essential properties of quantum chromodynamics (QCD) 
as the basic theory for the evalution of the equation of state for thermodynamical systems 
of elementary particles with strong interactions. The interest in the results discussed in this 
report has grown greatly due to the evaluations of the heavy ion experiments performed at 
CERN and BNL. Nevertheless, the details of our investigations are not directly dependent 
upon the outcomes of any given experiment. 

Any opening discussion of QCD calculations on the lattice begins with the original 
work of Kenneth Wilson [1] which is followed by the numerical simulations of Michael 
Creutz [2] as the important early developments in the field. In their early works the formu- 
lation of the lattice gauge theory and the suitable methods for evaluation were developed. 
Although we shall not discuss the actual numerical methods in detail, we shall try to indicate 
the approach. Since QCD is a gauge invariant theory, an important issue is the advancement 
of methods which uphold this property on discrete space-time points. Thereafter we can 
continue with a more general discussion of the methods and results [3], which we today can 
associate with the very extensive numerical evaluations in QCD at finite temperatures and 
densities. In the following sections of this report we shall discuss first the numerical evaluations 
for pure gauge theory at finite temperature with only heavy static quark sources, following 
which we include the dynamical quarks in the thermodynamics [4, 5] . Throughout this work 
we shall concentrate on writing the lattice QCD results in terms of the actual physical 
variables which are not simply the ratios of the measurable thermodynamical quantities. In 
this context a number of new quantities or variables have been introduced for the needs of 
these simulations on the lattice. In most cases these lattice quantities can be readily related 
to the corresponding quantites appearing in the usual quantum thermal field theory. In the 
cases of the thermodynamical quantities or densities it has been a custom to set in ratios 
of the densities for the sake of the computation instead of the actual physical variables of 
the continuum field theories. In this context we must explain why the defined [6, 7] lattice 
quantity A(T), which is commonly called the interaction measure written as (e — 3p)/T 4 . 
It replaces the actual physical quantities in the equation of state of the form e — 3p, where 
these are defined as the (internal) energy density e and the pressure p. From the lattice point 
of view the pressure ratio p/T A may be easily and accurately computed by direct integration 
techniques [8]. However, the quantity A(T) involves the renormalization group beta-function 
which was somewhat later for SU{2) gauge theory well computed [9]. This work led to a useful 
general procedure for the computation of the thermodynamical functions [10] for lattice gauge 
theories with SU(N C ) symmetries with different numbers of colored states, which of particular 
physical interest is clearly N c = 3. 

Above and beyond the pure numerical computations is fact that we must declare 
how these actual physical quantities may be related to the expected phase transitions which 
have been generally believed to take place. This consideration has been an active subject 
over the last quarter century [11] under the title of Quark Matter for which the disordered 
phase is oftentimes called the Quark Gluon Plasma (QGP). The relationship to the lattice 
computations as well as to perturbative QCD has been previously discussed [12]. In the very 
simplest picture an ideal gas of hadrons is converted into an ideal gas of quarks and gluons. 
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This picture fails to explain in any way why such a conversion should or even could take 
place. The first obvious improvement to this highly oversimplified model is the introduction 
of the bag constant B, which offers a means of holding the quarks and gluons together inside 
the hadrons (see Appendix A: Phenomenological Models for Confinement). The MIT bag 
model [13] provides a constant vacuum energy density e = + B and a vacuum pressure of 
p — — B for the hadrons in all directions. This special assumption makes the ground state of 
the hadrons more favorable than that of the free quark-gluon gas. Out of these two conditions 
on e and p the trace condition on the hadronic ground state average of the energy momentum 
tensor (Ofy = 4£> has been previously presented for the bag model [14]. We recall that the 
presence of a finite trace yields a breaking of the scale and conformal symmetry. Thus it is 
with these special properties of the bag model that the statistical theory for the hadrons with 
the inclusion of the strong interactions really starts. 

For QCD at low energies there arises another important issue which comes out of 
the spontaneous breaking of the chiral symmetry [15, 16]. It is known from atomic physics 
with many particles that the Nambu-Goldstone modes appear in the ground state when there 
is a spontaneous symmetry breaking 1 relating to a conserved current. In many known cases 
like the ferromagnet and the superconductor these modes remain with the quantum state even 
at finite temperatures. From the standard model of elementary particle physics the relation 
of this phenomenon to the lattice data will be discussed more extensively using a model in a 
later section. In this manner these modes are generally understood [16] using the sigma model 
with the SU (2) chiral symmetry for very low energies. Within small corrections to isospin 
symmetry the light quarks have the pion as the approximately massless Nambu-Goldstone 
mode, for which the vacuum expectation value (qq)o has the same value for both of the 
light quarks. If the masses of the light quarks were taken to be zero, the pion would be 
the true Goldstone boson with M n = 0. In the low energy limit with quarks of masses 
m u and the pion mass can be written as M n = (m u + m^)6 and (qq)o = —F^b , where 
bo is a positive constant of dimension energy and F% is the square of the pion decay constant [16]. 

Along this same line another important example of this type of approach was used 
long ago in the discussion of the evidence for the scalar meson dominance by Peter G. O. 
Freund and Yoichiro Nambu [18]. Out of this discussion we note the possibility of mesonic 
bound states at high temperatures. In this approach the trace of the energy momentum tensor 
is coupled through the Klein-Gordon wave equation to a single massive scalar field. In this early 
work they provide an effective Lagrangian for the dominance of the scalar meson. We shall 
develop this topic further in a later section. Another model with spontaneous symmetry break- 
ing which relates to the low energy properties of QCD is the Nambu- Jona-Lasinio Model [19,20]. 

In the course of this article we shall look at some of the work of the more recent 
past involving the transitions between the hadronic states and the the long speculated 
QGP. It involves both the restoration of the chiral symmetry as well as the deconfinement 
transition [21, 22]. We start our ideas with an development closely tied to some earlier 
work [14] on the hadron to quark-gluon phase transitions motivated by the relativistic 
heavy-ion collisions. Here we expect the properties of the quark condensate to have a 
significant part in the development of any new phase. The rules for the scaling properties 

1 A very recent discussion on this subject by Francois Englert with the title "Broken Symmetry and Yang- 
Mills Theory" can be found in Gcrardus 'tHooft's collection [17]. 
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then relate to the surrounding medium as proposed some time ago by Gerry Brown and 
Mannque Rho [23, 24] for the ratios of the decay constants and the masses in the presence 
of the surrounding medium to those of the vacuum. In this context we shall look into the 
use of the effective Lagrangian approach which will be used primarily in relation to the 
restoration of the chiral symmetry. Furthermore, we note that some more recent work on 
the nature of the chiral restoration transition [25] has been performed. A related discussion 
arises with the chiral bag model [21,26], where the action is constructed in such a way that 
it is invariant under global chiral rotation (see Appendix A for more detail). This model 
is an extension of the usual bag model which had ignored the properties of chiral symmetry [21]. 

Next we note the importance of some work on the use of the QCD sum rules at low 
temperatures [16,22]. This work was then related to some earlier numerical lattice simulations 
at finite temperatures [27] which involved the structures of the the electric and magnetic 
condensates separately. Along this line we should mention the important distinction between 
the "hard" and "soft" glue arising from the type of breaking of the different symmetries. 
Now it is quite necessary to note that there are two different types of symmetry breaking 
involved-that mentioned above as the spontaneous breaking involving the chiral symmetry and 
that which appears as the anomalous breaking of the conformal symmetry. The spontaneous 
breaking of chiral symmetry already appears in the hadronic ground state with the destruction 
of the chiral invariance in the axial current. It involves the operator average in the hadronic 
ground state of the form {ip q i(j q ) which has lower field dimension than the lagrangian density 
since the quark- ant iquark pair ip q ip q alone has the operator dimension three. However, the 
anomalous breaking of the scale and conformal symmetry arises from the square of the gluon 
field strength tensor, which we write symbolically as (G 2 ) in the hadronic ground state. 
It possesses the field dimension four. In the more general context it appears that with the 
loss of conformal symmetry, which we shall later see relates to the gluon condensate itself. 
It is never really restored even at very high temperatures. In the finite temperature field 
theory [4] another type of breaking occurs from the renormalization group equation at finite 
temperatures. The effect of the finite temperature renormalization first takes out the vacuum 
gluon condensate. Then, as it was clearly stated by Heinrich Leutwyler [28], it then continues 
to decondense with the increasing temperatures. This particular situation we shall discuss in 
the following paragraphs. 

Here we discuss further the ideas concerning the two different types of symmetry, 
which have been recently related to the problem of the mass in the mesonic bound states [29], 
which we shall discuss later in relation to the scalar meson dominance [18]. The study of the 
breaking of the chiral symmetry in gauge theories has had a very long history in quantum field 
theory. It had already arisen in other models well before QCD. The anomalous electromagnetic 
decay of the neutral pion, that is ir° — > 27, served as a major problem even before 1950. Its 
relationship to the study of anomalies was realized later in 1969 from spinor electrodynamics 
by Stephen Adler [30,31] and the work of John Bell and Roman Jackiw on the nonlinear sigma 
model [32]. It was later recognized as the anomalous breaking of chiral symmetry in QCD, 
which is usually known as the chiral or the ABJ anomaly [33-35]. It represents the anomalous 
divergence 2 of the axial current arising out of an axial Ward identity [34]. In contrast to 
the conformal or trace anomaly, which will be very essential to all further discussions of the 

2 A very interesting discussion of both the chiral and the trace anomalies has been recently written by Stephen 
Adler entitled "Anomalies to all orders", in reference [17]. 
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physical equation of state in QCD, the chiral or the ABJ anomaly had always a higher rating 
in its acceptance because of the clear advantage from the already well known experimental 
verifications using the decays of the 7r° as well as that of such hadronic processes as the K + 
going into pions [34]. 

The discovery of anomalous terms appearing as a finite value of the trace of the 
energy momentum tensor was pointed out as a result of nonperturbative evaluations in 
low-energy theorems [42] many years ago. Furthermore, it was also somewhat later realized 
how this factor arose with the process of renormalization in quantum field theory which 
became known as the trace anomaly [44,46,47] since it was found in relation to an anomalous 
trace of the energy momentum tensor. The basic idea of the relationship between the trace 
of the energy momentum tensor and the gluon condensate has already been studied for 
finite temperature by Leutwyler [28] in relation to the problems of deconfinement and chiral 
symmetry. The starting point of this work begins with a detailed discussion of the trace 
anomaly based on the interaction between the Goldstone bosons in chiral perturbation theory. 
Quite central to his investigation is the role of the energy momentum tensor averaged over 
all the states, whose trace is directly related to the averaged gluon field strength squared. 
Here it is important to state that the averaged total energy momentum tensor T^ U (T) can be 
separated into the vacuum or confined part, 6q V , involving only the temperature independent 
states in the average, and the finite temperature contribution 6^ U (T) as follows: 

T""(T) = 9^ + 9^ u (T). (1) 

The temperature independent part, 6q V , has the standard problems with infinities of any ground 
state, which has been previously discussed [36] in relation to the nonperturbative effects in QCD 
and the operator product expansion. In the following analysis we shall start our discussion 
with a bag type of model [13, 14] as a means of stepping around these difficulties with the 
QCD vacuum since at this time we are primarily interested in the thermal properties 3 . The 
finite temperature part, which clearly vanishes at zero temperature, has no such problems 
with the divergences. We shall discuss in the following sections of this report how at finite 
temperatures the diagonal elements of 6^ V (T) are calculated in a straightforward way on the 
lattice. Furthermore, the trace 0^(T) is connected to the thermodynamical contribution to the 
internal energy density e(T) and pressure p(T) for relativistic fields as well as for relativistic 
hydrodynamics in the following simple form: 

0£CO= e (T)-3p(T). (2) 

The gluon field strength tensor including the coupling g is denoted by G^ u , where a is the color 
index for SU(N C ). Thus the basic equation for the relationship between the gluon condensate 
and the trace of the energy momentum tensor at finite temperature was written down by 
Leutwyler [28] using the trace anomaly. Leutwyler's equation takes the following form: 

(G 2 ) T = (G 2 ) - e^Tl (3) 

for which the brackets with the subscript T mean thermal average and the zero the ground or 
confined state average. The renormalized gluon field strength tensor is squared inside of the 

3 For a discussion of the process of gluon condensation in relation to confinement see the contribution of 
David Pottingcr in the collection on the statistical mechanics of quarks and hadrons [11]. 
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brackets which is then summed over all the colors to yield 



G 2 = ~^G^G%. (4) 



The renormalization group beta function (3(g) in terms of the coupling may be generally written 



MS 



M =V^f= -7^( lliVc - 2N f^ 3 + 0(g 5 )- (5) 

Out of these relationships Leutwyler [28] has calculated the trace of the energy momentum 
tensor at finite temperature for two massless quarks using the low temperature chiral expansion . 

The most immediate generalization of Leutwyler's equation Q has been previously 
considered in an earlier work [64]. In the presence of massive quarks the averaged trace of the 
energy-momentum tensor takes the following form from the trace anomaly: 

^ = rn q $ q ip q ) + (G 2 ), (6) 

where m q is the renormalized quark mass and ip q , ip q represent the quark and antiquark fields re- 
spectively. As an operator relation this equation (unaveraged) would not make any sense, since 
these three operators have different operator dimensions and also carry different symmetries 
therein. We include with these averages the renormalization group functions (3(g) and ~f(g,m), 
which appear in this trace from both the coupling and mass renormalization processes. This 
averaged form holds for both the confined and temperature dependent structures. We shall 
use these properties in relation to the massive dynamical quark lattice simulations in Section III. 



We start with a discussion of previous evaluations [37] from an earlier collaboration 
with Graham Boyd. Originally the motivation for this work was just to study the pure 
lattice gauge theory data for both the SU (2) [9] and the SU (3) [10] simulations, which had 
at that time been recently finished in Bielefeld. Of particular interest at that time was 
the relation of the equation of state to the pure gluon condensate above the deconfinement 
temperature, which had been in both cases very accurately computed and denoted as T c . As 
it had been stated previously by Leutwyler [28], we did, indeed, find that with increasing 
temperatures the pure gluon condensate was unbounded from below over its range of negative 
values . This observation was quite contrary to the then commonly accepted ideal gas models 
which were supposed to appear at high temperatures because of "asymptotic freedom" in 
QCD. Actually in our present understanding we know that the important property of QCD, 
asymptotic freedom, is already in the renormalization group beta function whose stable 
ultraviolet fixed point determines the critical behavior at T c . During the time that this earlier 
work [37] was being carried out there appeared some numerical lattice computations from the 
MILC collaboration [48, 49] which included two light dynamical quarks. At the same time in 
Bielefeld [51] there were similar computations for four flavored not so light dynamical quarks. 
These results we included at the end of our work [37] (see our Figure 4 where these results are 
compared to a rescaled SU(3) curve). We noted, but could not explain why, that the heavier 
Bielefeld data followed the rescaled pure gauge curve, while the MILC data [48,49] remained 
well above it. There was then, as there is now [50], the problem that this data for the two light 
quarks did not go to a high enough value in the temperature range to make such comparisons 
clear. Although this particular work [37] never got published, it did, nevertheless, serve as a 
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starting point in further works [52,64,65]. 

An important result of the investigation of numerical lattice simulations was already 
stated by Leutwyler [28] with some of the previously cited numerical work. This was the fact 
that the dilatation current, given by = x u T^ u , and the four conformal currents K afl (x) are 
not conserved quantities even at very high temperatures. The four conformal currents [35] are 
given by 

= (2x a x u - g«x 2 )T» v . (7) 

The form of the equations for these currents was first written down by Erich Bessel-Hagen [66] 
as the result of a seminar series in 1920 at Gottingen led by Felix Klein . The equation for the 
dilatation current takes the form 

d,D» = T£. (8) 
A similar equation can be written for the divergence of the special conformal currents 

d^K^ = 2x a T^, (9) 

It is important to note that the nonconservative aspect of both these currents D^(x) and 
K fia (x) relate directly to the fact that remains finite. We shall discuss some aspects of 
these currents with special solutions as well as the corresponding differential and integral forms 
in Appendix B. However, in this report it is our main interest to discuss these quantities from 
the point of view of lattice data. 

We shall show in the very next section how the actual computed thermodynamical 
functions appear as a function of the temperature in physical units. Thereby in the next part 
of this report we discuss how to use the lattice data for pure gauge theories 4 . In particular 
the method for the evaluation [6, 7] of the lattice quantities as the ratios p/T 4 , e/T 4 and 
A(T) is the real starting point for the discussion of the thermodynamical quantities leading to 
the form of the equation of state. It symbolizes the change of scale breaking from the pure 
vacuum contributions to the dominance of the high temperature region in the thermodynamics. 
Furthermore, it is significant that just above the deconfinement temperature the relative 
breaking is the largest as seen from the peak in A(T) with the gradual decline thereabove. 
The actual equation of state e(T) — 3p(T) shown in the next section demonstrates this fact 
much less dramatically just above Td since its rise is very sharp there which at much higher 
temperatures then simply slows down. In the following section we consider the inclusion of 
dynamical fermions in the numerical computations. Although we shall generally mention some 
of the studies with finite chemical potentials on the lattice as well as later consider some 
special cases as examples, we will not take up the details of the investigations here. Instead, in 
the next section we shall provide numerical analyses of various cases for the finite temperature 
behavior of the quark and gluon condensates relating to the breaking of chiral and conformal 
symmetries. The last major section attemps to put together these numerical results with the 
contemporary understanding in nuclear and high energy physics at high temperatures. A short 
section with concluding remarks is followed by two appendices with added details. 



4 In this report we shall only introduce the notation needed for the lattice evaluations. Here we will not 
discuss the usual field theoretical notation for QCD. We shall follow the standard texts on gauge fields such 
as [16,20,56,57]. For discussions of the ideas and their development see [17]. 



8 



II. Thermodynamics of Gauge Theory on the Lattice 



Here we begin by defining the physical quantities in terms of the lattice varibles [f] 
used in the following parts for the lattice gauge computations. First we discuss the ther- 
modynamics of the pure Yang-Mills fields as it is computed on the lattice [2, 3] in the 
canonical ensemble in the formalism necessary for finite lattice sizes. A basic fact of the pure 
gauge theory is the existence of a transition temperature often known as the deconfinement 
temperature or For the SU (2) invariant gauge theory this transition temperature is often 
called the critical temperature T c since this transition is of second order [6,9]. However, 
for the SU (3) invariant gauge theory it is a first order [7, 8, 10] phase transition which goes 
between the "confined" and the "deconfined" states. In this case the critical temperature 
acts as the highest temperature where there is a distinction between these phases 5 . In 
this part we go more specifically into the numerical results of the actual lattice computa- 
tions for the thermodynamical functions using lattice gauge theory at finite temperatures [9,10]. 

II. 1 Lattice Thermodynamics 

As it is usually done in statistical physics, we start with the canonical partition func- 
tion Z(T, V) for a given temperature T and spatial volume V. From this quantity we may 
define the pressure p for large homogeneous systems in thermal equilibrium through its relation 
to the free energy density / as follows: 

V = ~f = ^\nZ(T,V). (10) 

The volume V is determined by the lattice size N a a, where a is the lattice spacing and N a 
is the number of steps in the given spatial direction. The inverse of the temperature T is 
determined by N T a, whereby N T is the number of steps in the (imaginary) temporal direction. 
Thus the simulation [8-10] is done in a four dimensional Euclidean space with the given lattice 
sizes iVj x N T , which gives the volume V as (N a a) 3 and the inverse temperature T' 1 as N T a 
for the four dimensional Euclidean volume. 

In the early lattice simulations for the SU(2) gauge theory one took for the actual 
lattice sizes up to N a = 24 in the spatial directions and the temporal sizes iV T = 4, 6 [7,38]. 
Later for the SU(3) lattice simulations the spatial values were taken as N a = 16, 32 with the 
temporal sizes N T = 4, 6, 8 [10], which we shall show here. In general for SU(N C ) gauge theory 
the lattice spacing a is a function of the bare gauge coupling (3 defined by 2N c /g 2 , where g is 
the bare SU(N C ) coupling. Thereby this function fixes both the temperature and the volume 
at a given coupling. Now we write P a<T as the expectation value of, respectively, space-space 
and space-time plaquettes in terms of the link variables U 

P*,r = 1 - ^-ReiTriU^uM)) (H) 

for the usual Wilson action [10]. These plaquettes P CT)T may be generalized to the improved 
actions on anisotropic lattices [39] for SU(2) and SU(3). For the Wilson action we define the 

5 The presence of a transition temperature to a high temperature phase was suggested many years ago from 
renormalization group arguments. This phase was regarded as where the gauge coupling became weaker [67] 
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parts 5*o = QPo on the symmetric lattice and St = 3(P CT + P T ) on the asymmetric 
lattice N^xN T . We now proceed to compute the free energy density ratio defined above in the 
equation (fTUj) by integrating these expectation values as 

P 

J[ ' !) - -N* [d(3'{S - St], (12) 



ji4 

h 

where the lower bound (3q relates to the constant of normalization. At this point we should 
add that the free energy density is a fundamental thermodynamical quantity from which all 
other thermodynamical quantities can be gotten. Also it is very important in relation to the 
phase structure of the system in that the determination of the transitions for their order and 
critical properties as well as the stability of the individual phases are best studied. The integral 
method [8] for the computation of the pressure ratio yields as a result from the equations (fTUj) 
and (J12)) given by 

P 

Yl = |l + ^Jd(3'[S - S T ]. (13) 

Po 

The integral method provides an approach for the evaluation of the pressure ratio which is free 
from many of the problems arising in the evaluation of the ratios of other thermodynamical 
functions as well as assuring that the value of the pressure ratio always be positive as 
distinguished from many of the earlier evaluations. For a general discussion see reference [3]. 

Next we define the lattice beta function in terms of the lattice spacing a and the 
coupling g as follows: 

da- 2 

(3(g) = -2iV c a-^. (14) 
The dimensionless interaction measure A(T) discussed in the Introduction [6] is then given by 

A(T) = iV 4 /%) [S - S T ]. (15) 

The crucial part of the more recent lattice gauge calculations is the use of the full lattice beta 
function, (3(g) in obtaining the lattice spacing a, or the scale of the simulation, from the coupling 
g 2 . Without this accurate information on the temperature scale in lattice units it would not 
be possible to make any claims about the behavior of the gluon condensate discussed in detail 
in Part IV. The interaction measure A(T) is the thermal ensemble expectation value given by 
(e — 3p)/T 4 . Thus because of equation (J2J) above the trace of the temperature dependent part 
of the energy momentum tensor, here denoted by 9fi(T) is equal to the expectation value of 
A(T) multiplied by a factor of T 4 . This physical quantity may be directly computed [37,52] as 
a function of the temperature as 

0£(T) = A(T)-T 4 . (16) 

There are no other contributions to this trace for the pure gauge fields on the lattice. The 
heat conductivity is zero. Since there are no other finite conserved quantum numbers and, 
as well, no velocity gradient in the lattice computations, hence no contributions from the 
viscosity terms appear. For a scale invariant system, such as a gas of free massless particles, 
then the trace of the energy momentum tensor in equation (|16|) is clearly zero. However, any 
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system that is scale variant, perhaps from a particle mass, has a finite trace, whereby the value 
of the trace then measures the magnitude of scale breaking. These results are shown in figure ^ 



We now discuss a few other important properties of the pure gauge theories on the 
lattice. The order parameter is often taken as the lattice average of the Polyakov loop [3] , 
which in many cases turns out to give more information as the order parameter than does the 
general Wilson loop. The defining operator for the Polyakov loop is given by 



1 v T 

,1 J -X4 = l 



(17) 



where the index 4 stands for the euclidean time direction. Then the actual Polyakov loop is 
defined by the expectation value as 



Since a vanishing expectation value is generally a signal of an exact symmetry, one oftentimes 
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Figure 1: We compare the computed values of trace of the energy momentum tensor in the 
units [GeV 4 ] for the lattice gauge theories £77(2) and SU(3) as indicated. The values for the 
deconfinement temperature are Ta = 0.290,0.264 GeV for SU{2) and SU{2>), respectively. 
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defines the absolute value of the quantity L as 

z = ^<iE/®i>, ( 19 ) 

which serves as an approximate order parameter. For the pure gauge theories one can determine 
the critical couplings (3 C from the analysis of the Polyakov loop susceptibility xl, 

Xl = K[(L'} - (L) 2 ] , (20) 

which we shall later discuss in more detail in conjunction with the dynamical quarks in the next 
section. From this calculation one can accurately determine the value of the deconfmement 
temperature T d . 

II. 2 Thermal Field Theoretical Evaluation 

In a very recent work Zwanziger [53] has carried out an analytical determination of the 
properties of the equation of state for the pure gluon plasma at high temperatures well above 
the deconfmement temperatures. He uses the Gribov dispersion relation [54] as a means of 
suppressing the infrared modes when gauge equivalence is imposed at the nonperturbative 
level, which takes the form 

E(k) = ^k 2 + M A /k 2 . (21) 

After a brief analysis of the trace anomaly at finite temperature [28] as well as comparison 
with the SU(3) lattice data [10] Zwanziger derives a form of the trace of the energy momentum 
tensor 

0£(T) = LT + A, (22) 

where A is a temperature independent term of order 0(1) in the units [Gel/ 4 ] and L is an 
integral of the form [53] 

POO 

L = 3(N 2 C - 1)/tt 2 / dkk 2 \n[E{k)/kl (23) 
Jo 

This integral can be exactly evaluated for the Gribov dispersion relation (|21j) yielding the 
following result: 

L = (N 2 - 1)M 3 /ttV2. (24) 

Although this result with a linear growth in the temperature for 9£(T) had been already 
observed [37, 52] in the lattice data, it had not been previously derived analytically for this 
energy spectrum [53]. If one looks carefully at the figure [U for 9jt(T), one sees well above 
the deconfmement temperatures Td that both the curves begin to slow up to an almost 
linear behavior as a function of temperature T even for the temperatures of around 2T^. 
Nevertheless, there still are some quite gradual changes of curvature upwards until almost 
3Td or about T = 0.8GeV, whereupon, it appears for higher temperatures to be quite 
linear [37]. Furthermore, one sees a new scale M in the equation (|21|) which has the units 
[GreV]. According to Zwanziger [53] it arises from the infrared regulator, which can be 
estimated to be around 0.7GeV^. From a more general standpoint what is very interesting 
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about Zwanziger's analytical calculation of Ofi(T) is that the linear result for the temperature 
arises from the regulation of the infrared behavior of QCD which dominates over this high 
temperature domain. It is by now very well known that there are serious infrared problems 
in the high temperature perturbation theory [55], which are not easily dealt with using the 
standard perturbation theory. These problems indicate that the asymptotic freedom may 
only be strictly valid in the region of short distances and times [3]. We already know that 
for QCD the asymptotic freedom has appeared in the beta function as a property of the 
ultraviolet structure. This property of QCD provides a stable ultraviolet fixed point in the 
beta function in contrast to quantum electrodynamics (QED) which has a stable infrared 
fixed point [20,56,57]. However, this fact means that QED is infrared stable in contrast to QCD. 

II. 3 Numerical Evaluations of Physical Quantities 

The numerical evaluation of the equation of state at finite temperature for strongly in- 
teracting quarks and gluons has long been the main objective for lattice simulations [2,3,5]. 
The computed pressure ratio p(T)/T 4 for the pure SU(3) lattice gauge theory [10] is shown as 
a function of temperature in their Figure 4 with three different values of N T in the Euclidean 
time direction. They remark concerning the sizes of the lattices for the different cases N a 3 xN T . 
In what follows we take the spatial sizes N a a and temporal sizes N T a whereby a is the lattice 
spacing between points on the lattice. In these particular computations [10] the spatial step 
numbers have the two values N a = 16, 32 while the temporal number have the three values 
N T = 4,6,8, each of which has been explicitly studied [10]. We should note that in their 
Figure 4 the value for the "critical" deconfinement temperature called T c for the different 
lattice sizes has been been evaluated from the renormalization group beta function on the 
lattice. The attained value for T c from an extrapolation to the continuum limit is given by 
T c / \fo = 0.629(3). The string tension is taken for their simulations as = 420MeV, which 
results in a "critical" temperature of about T c = 2QAMeV. 

The actual simulations [8] for the pressure ratio in SU(N C ) around T c was proposed 
by using the integral method, which prevented the undesirable effects in some earlier 
evaluations of the pressure, in which the pressure ratio appeared to be negative near the 
critical coupling [58, 59]. Furthermore, the computational effect of the lattice anisotropy 
can be very accurately computed by using the integral and differential methods for the 
anisotropy coefficients leading to a much higher resolution. This program has been discussed 
in detail more recently [39] by showing how these lattice computations for both SU(2) and 
SU(3) are carried out. Similarly the energy density e(T) is calculated from the pressure p 
and the interaction measure A(T) to arrive at the form A(T)-T 4 + 3p as a function of 
temperature. This method was briefly mentioned earlier. These numerical results are shown 
in a similar plot in their Figure 6 with the same basic parameters as those for the pressure 
ratio [10]. In their Figure 6 we can see that the energy density ratio e(T)/T 4 as a function 
of the temperature T for pure SU (3) gauge theory rises much faster around the critical 
temperature T c than does the pressure ratio p(T)/T 4 in the previously mentioned Figure 4. 
It is evaluated from the interaction measure as a function of the coupling. In the sense of 
the thermodynamics the internal energy is generally used to thermally describe the state of 
the system. Thus we may well expect that its derived form as a density, e(T), should be 
more sensitive to the change of phase at T c . From a casual view of the number scale we 
are able to see that the energy density ratio curves in each case for N T always lie consider- 
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ably above that of the pressure- even above three times the corresponding numerical values [10]. 

Now we look explicitly at the equation of state of the pure gauge theory which is 
gotten directly from A(T)-T 4 . Thus this quantity which we have discussed extensively in 
the Introduction is simply related to the trace of the energy momentum tensor 9£(T) given 
above in the equation (J2J). Here we plot 6 the equation of state e — 3p as a function of T/T c . 
Furthermore, we want to point out carefully the continual growth of the equation of state for 
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Figure 2: The equation of state e — 3p for pure SU(3) gauge theory with the three lattice sizes 
in the physical units [GeV ]. The temperature is given as T/T c . 

pure SU (3) gauge theory as shown in the Figure 121 This fact arising from these numerical 
simulations shows a behavior that is quite contrary to many of the common speculations on 
the equation of state for the quark-gluon plasma. We can see here no obvious signs from any 
of these computations that the dependence of this quantity e(T) — 3p{T) decreases to zero at 
any temperature above the "critical" deconfmement temperature T c . In fact, one can reaffirm 
here with only very minor variations for the sizes of the different finite lattices that the high 
temperature linear T dependence theoretically calculated [53] is still quite well upheld. Thus 
it is safe to conclude that the pure SU(3) gauge theory in the computed range of temperatures 
remains a strongly interacting system of gluons- not an ideal ultrarelativistic gas! 

6 The author thanks Jiirgen Engels for pointing out errors in earlier graphs and replacing these plots in a 
corrected form from the original data [10]. 
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II. 4 Comparison of the Physical Quantities 



In this paragraph the results of all these computations are summarized in the last fig- 
ure in this part. In this Figure El we make an actual comparison using the physical units which 
are written below the curves in the square brackets. Here we can clearly see the differences 
between these computed thermodynamical quantities as functions [10] of the temperature 
T[MeV]. In addition to the above analysed quantities we have included the entropy density 
s[fm~ 3 } which is obtained from (e + p)/T, all of which are known. However, the entropy 
density s[fm~ 3 } shown as the upper curve on the right cannot be directly compared with the 
others because of the units. The equation of state e — 3p[GeV/ fm 3 ] as the lower curve on 
the right can be compared with the energy density e and the pressure p, since all are in the 
same units. One can see that for these physical quantities the relative growth of each at the 
deconfinement temperature T^. Above each one flattens out at different rates. This effect 
we can clearly see in the comparison for pure gauge theory in Figure E3 The growth of the 
equation of state is obviously very much slower than all the other thermodynamical quantities. 
This difference amounts to the comparison of a linear increase in the temperature to those 
with cubic or quartic powers when taken on a logarithmic scale. 



1000 r 




200 300 400 500 600 700 

T [MeV] 

Figure 3: The different thermodynamical functions of the temperature for pure gauge theory. 
The vertical line is the deconfinement temperature. Here we take for these comparisons the 
more useful density units rather than the pure energy units as in the previous figures. 
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There has been some recent high temperture work in perturbation theory [61] to higher 
orders in the coupling g, for which the pressure ratio to the Stefan-Boltzmann limit has been 
compared to the lattice results for the pure gauge theory [10]. Furthermore, a comparison 
between the perturbative results [62] up to the order O (g 6 ) and these lattice gauge simulations 
in Figure 121 for e — 3p shows a very good agreement above about 800MeV\ This perturbative 
evaluation 7 for e — 3p continues to increase well beyond the range of the lattice data in roughly 
the same way as the simulations [10]. 

II. 5 Discussion of Conformal Symmetry 

Specific investigations of the effective measure of conformal symmetry have been re- 
cently carried out on the lattice at high temperatures for pure SU(3) gauge theory. As we 
have already discussed above in the Introduction, the presence of a nonzero trace for the 
energy momentum tensor relates with the breaking of scale and conformal invar iance. We shall 
discuss more thoroughly in the Appendix B the mathematical nature of the related physical 
currents D^(x) and K^ a {x) which are relate directly to the fact that the trace of the energy 
momentum tensor remains finite. 

Then the question clearly arises concerning how close the given thermodynamical 
system is to achieving the scale and conformal invariances. A measure of the deviation 
from conformality is given by the expression C = (e — 3p)/e, which has been proposed 
by Rajiv Gavai, Sourendu Gupta and Swagato Mukherjee [60]. Clearly if the value of C 
were identically zero, then the conformal invariance of the theory would be fully upheld. 
However, they clearly find from their simulations [60] at both the temperatures of 2T^ and 
3Trf the regions with finite (nonzero) values for C. Thereby they show in a plot of the 
pressure ratio p/T 4 against the energy density ratio e/T 4 that all the numbers lie significantly 
below the line of conformality in the pressure where C = 0, which corresponds to the ideal 
ultrarelativistic gas. Then from their work we are able to see that all of the values of p/T 4 
and e/T 4 in their simulations are visibly removed from their ideal gas values for a lattice 
with those given quantum numbers of the spin and color. Thus for pure gluon system we 
would expect to find the Stefan-Boltzmann limiting values for these ratios to be 87r 2 /45 and 
87r 2 /15, respectively. Furthermore, the higher temperature points at 3T^ are distributed in 
a cluster which appears to be further from the line of conformality C = than were those at 2T^. 

Finally we conclude this part on the thermodynamics of the pure lattice gauge the- 
ory. In summary, for the gluon gas with strong interactions simulated for both SU(2) and 
SU (3) symmetries all the thermodynamical functions in physical units including the equation 
of state grow monotonically in the temperature. This clear statement from the computed 
numerical results is obviously quite contrary to the usual expectations from the ideal gas 
oriented theories. In some earlier investigations and collaboration with Graham Boyd [37, 52] 
we have shown further properties of these thermodynamical functions to possess a steady 
growth in the presently considered range of temperatures continuing to, at least, 1.5 [GeV]. 



The author thanks Mikko Laine for showing to him these numerical results from finite temperature 
perturbation theory. 
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III. Dynamical Quarks at finite Temperature 



In this part of the report we look at the full lattice QCD including the thermodynami- 
cal contributions to the equation of state arising from the thermal properties of the dynamical 
quarks. The presence of these quarks with multiple flavors changes radically the evaluation of 
the thermodynamical quantities. The changes are largely due to the presence of the broken 
chiral symmetry in the hadronic ground state of the colored quark fields. Furthermore, the 
restoration of the chiral symmetry at finite temperatures radically restructures the high 
temperature quark-gluon phases. The critical temperature T c for the chiral restoration is 
considerably lower than the deconfinement temperature Ta for the pure SU(N C ) lattice gauge 
theories discussed in the last section. The order parameter can be defined similarly to the 
Polyakov loop L in equation (j!8|) for the pure gauge theory except that one takes only the real 
part 3?e L when the dynamical quarks are present [68,69]. 

Before we consider the actual lattice data, we mention the thermodynamics of the 
ideal relativistic gas of rif flavored quarks with a given rest mass m at a temperature T. The 
quantities like the pressure and the energy density are well known, which gives the equation of 
state 

e(T) - 3p(T) = gm^TK-iim/T), (25) 

where Ki(m/T) is the modified Bessel function of the second kind. The statistical degeneracy 
factor g has the value 3n//7r 2 . For a fixed quark mass m and at temperatures large compared 
to the mass energy then the right side of (J25j) becomes just gm 2 T 2 . Then the equation of state 
of a pure massive gas grows quadratically in both the temperature and the mass in the high 
temperature limit. 

III.l Equation of State for two light Quark Flavors 

First we consider the equation of state for two light dynamical quarks which appears as 
a plot similar to the figure 121 for the pure SU(N C ) gauge theories in terms of the trace of 
the energy momentum tensor for the equation of state 9fi(T) = e(T) — 3p. In this new 
figure for the MILC97 data [50] we have shown the equation of state using equation (|16|) in 
terms of the physical units [GeV/ fm 3 ]. Furthermore, in this representation we are able to 
compare our results directly with those in their Figure 7 for the MILC data [50] which shows 
the interaction measure written as (e — 3p)/T 4 against the QCD lattice coupling 6/g 2 over the 
range of values from 5.36 to 5.50 with the two mass values am q = 0.0125 and am q = 0.0250, 
which are expressed in terms of the lattice spacing. 

It is important that we note here for the case of the light dynamical fermions that 
even though the pressure ratio is computed in a similar way to the pure lattice gauge theory 
starting with the equation (|10|). it, nevertheless, has the predominant effect of the average of 
the product of the antiquark ^ q and quark ip q fields written as (ip q ip q ). Thus the equation for 
the pressure ratio (|13j) gets modified to contain these contributions. Since we do not present 
here the numerical values for the pressure, one should go to the original MILC97 data [50] 
for the exact form of the pressure equation, which, however, uses slightly different notation 
in the integration method as applied to the free energy ratio. For this reason the interaction 
measure has an additional contribution beyond just the plaquette terms of the pure lattice 
gauge theory as given in equation (|15|). This new term arises from the mass renormalization 
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on the antiquark-quark terms (ip q tp q ). These two contributions together give the interaction 
measure A m (T) with the the quarks of mass m q 

A m (T) =n% 0{g) [3 - S T ] + ^am q ) [$M - (26) 

The beta function (3(g) is given similarly to that in equation (|14j) .However, the renormalization 
group gamma function for QCD [56] can be defined on the lattice 7(am 9 ) as 

d(am a ) , . 

The presence of both these terms in A m (T) carries very important consequences in the 
equation of state as can be seen in the figure EJ The behavior between around lOOMeV up 
to just below 250MeV for the MILC97 data [50], that is near to and just above T c shows a 
decisively different change in the equation of state from that of the pure SU(3) gauge theory 
just above as seen in figure El Although the pure gauge theory shows some variation in 
its rise above T d depending upon the lattice sizes [10], after about 300MeV its increase is 
reduced to an approximately linear rate [53]. In the later paragraphs we shall contrast both 
these results with some newer lattice computations for different numbers and masses of quarks. 
First, however, we shall discuss a little more carefully the structure of chiral symmetry breaking. 
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III. 2 Chiral Symmetry and Dynamical Quarks 



In the presence of dynamical quarks another symmetry becomes important- the chiral 
symmetry. When the quarks have masses, this symmetry is automatically broken. The chiral 
symmetry is a property of the two different representations of SL(2, C) denoted by 2 and 2* 
arising for the Dirac spinors in the Weyl representation [20]. It is the presence of the quarks' 
mass terms in the Dirac equation that formally breaks the chiral symmetry. This comes 
formally out of the nonconservation of the axial current j£ as discussed [20,35] relating to the 
triangle diagrams, such that the chiral anomaly for QCD takes the form 

d$ = k + ^~ 2 G^G; u , (28) 

where k\ is a constant and j$ is a pseudoscalar contribution. This situation has important 
implications in the case for finite temperatures where for T sufficiently high the chiral symme- 
try is restored in the small mass or chiral limit, m q — ► 0. We shall discuss the implications 
of this both from the theoretical side and the numerical side where a finite small mass is present. 

We now look at the chiral condensate at finite temperatures using chiral perturba- 
tion theory. The low temperature expansion for two massless quarks can be written [28] in the 
following form: 

' ' ~ 1 - o™ ~ - ^™ ^ + ^(T 8 ) + 0(ex P — ), (29) 



where F n is the above mentioned pion decay constant and the scale A q is taken as approxa- 
mately 0A70GeV. Leutwyler has shown that this expansion up to three loops remains very 
good at least up to around O.lOOGeV. At very low temperatures the probability of finding 
any given excited mass state is related to the exponentially small correction, which in this 
case has a very small value. As the temperature becomes higher, the number of different 
states begins to grow 8 . Nevertheless, at sufficiently low temperatures the excited states 
may be regarded as a dilute gas of free particles since the chiral symmetry supresses the 
interactions by a power of T of this gas of excited states with the primary pionic component [28]. 

Upon approaching the chiral symmetry restoration temperature T x the picture changes 
drastically. At this point the ratio T/F^ is considerably greater than unity. It is here where 
one expects the chiral condensate to be very small or to have totally to have vanished. This 
effect has been studied recently numerically [40] for two light flavors at finite temperature on 
the lattice. The results of this simulation is shown for ('0gV , ?)o> which we simply write 

as the quark condensate ratio, (ipip), in the following two plots for the restoration of chiral 
symmetry. We show this quark condensate ratio as a function of the coupling (5 for the range 
where the chiral symmetry is mostly restored [40] . The Figure |5] shows this ratio for two light 
quarks with a mass in lattice units of 0.02 on a lattice of size 16 3 x 4. We remark that at 
the value for the coupling 5.24 the chiral symmetry for the two light quarks with the lattice 
mass of 0.02 is already about 20% restored. While at the upper value of 5.33 it is still only 
about 80% restored. Thus we see that the chiral limit has not in this case been reached. The 



8 We note that here the behavior of the quark condensate (t/'gV'g) is no longer dominated by the low energy 
meson states like the pions and kaons. As the energy increases the number states begins to grow exponentially 
as would be indicated by the Hagedorn spectrum [87], which leads to an expected problem with this type of 
series at high temperatures. We will mention this situation later in the report. 
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Figure El comparess how the different mass values shown from left to right of 0.02, 0.0375 and 

0. 075 depend upon the given lattice sizes, which are 8 3 x 4, 12 3 x 4 and 16 3 x 4. We should 
also notice how the larger mass values slow the restoration down, which corresponds to moving 
the transition T x to higher temperatures or even eliminating it altogether as indicated by the 
flatness of the curves. 

III. 3 Thermodynamics with Dynamical Quarks 

The main quantities which were analyzed here were the various susceptibilities: 

1. The Polyakov loop susceptibility we have defined earlier in the last section in the equation 
fl2()(l . We will use now in connection with two other forms of the susceptibilities which are to 
follow: 

2. The magnetic or chiral susceptibility; 

x™ = v^l^ lnZ{T > n (30) 




5.24 5.26 5.28 5.3 5.32 5.34 

Figure 5: We show < ipip > as a function of the coupling (3 for the quark mass in lattice units 
ma= 0.02, which is normalized to the vacuum value of the chiral condensate. 
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3. The thermal susceptibility; 



Xe 



lnZ(T,V). 



(31) 



One compares the critical properties of xl, Xm and xe in order to establish the value of T x 
and its critical properties in the chiral limit where m q — > 0. For the moment we use T x for 
the chiral restoration temperature in contrast to the critical temperature T c for the dynamical 
quark simulations. However, in numerical simulations TOj must be taken to be finite- this means 
that one must use various different small values of on the different sized lattices N% x N T . 
This procedure uses the lattice data to find the values around the peak of the susceptibility Xm 
at T x for the smallest masses, with which one can determine the critical structure. A careful 
determination of the topological susceptibility relating to the chiral current correlations can be 
related to the square of the topological charge Q\ in the chiral limit [41], such that 



n 



m 



(Q 2 T ) = V(W) 



m— >0' 



(32) 



where n/ is the number of flavors. Thus from these susceptibilities one can arrive at the quark 
condensate ratio {fiif)). However, in this computation it is a major problem to properly set 
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Figure 6: We show the same physical quantity with the different values of the quark masses 
0.02, 0.0375 and 0.075 from left to right [40]. 
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the temperature scale for small lattices with finite masses. The plots in the figure El and El are 
made with the coupling (3 which may be compared with pure £77(3) on one side and the two 
flavor dynamical quark simulations on the other [50]. In the case of pure £77(3) the critical 
coupling (3 C for a 16 3 x 4 lattice has the value [10] of about 5.70, which is considerably larger 
than the values of (3 shown in the figure Eland El However, for the two light flavored dynamical 
quarks [50] the value of (3 C is around 5.40, which is still somewhat above these values shown in 
the two figures. 

Here we have investigated the properties of the chiral symmetry restoration for the 
quark condensate (ipqipq) alone at various values of the coupling (3. We have also shown in the 
figure El how the coupling shifts to higher numbers for larger values of the mass. These results 
can be compared to the MILC97 data [50] with much lighter quark masses which, nevertheless, 
stays mostly in the same range of the couplings (see in the reference [50] their figure 4.) We 
can see there that this data is extended into a higher range of coupling to get considerably 
smaller values of the quark condensate. These values we shall use in the next section for the 
mass contribution of the quark condensate at higher values of the temperature. We can then 
compare this effect in the chiral limit. However, here it is very difficult to immediately go over 
to a physical temperature scale in the same way as in the previous section for the pure gauge 
or gluon system. In what follows we shall look into the gluon condensate in the presence of 
dynamical quarks. Here we know that the presence of the quark masses are an immediate 
cause of scale symmetry breaking which of course change the scale of the system. This in turn 
changes the beta function as well as adds a term due to the mass renormalization. Thus the 
renormalization group equations are changed accordingly. This effect we shall discuss more 
thoroughly in the following. 

III. 4 Equation of State with Different Quark Flavors 

Next we look at the thermodynamical functions including the pressure p{T), energy 
density e(T) and the equation of state e(T) — 3p(T) all in the physical units [GeV/ fm 3 } for 
the case of somewhat heavier dynamical quarks with different numbers and types of flavors. 
These different cases have been worked out in the doctoral thesis of Andreas Peikert [72] 
at the Universitat Bielefeld for the lattices of sizes 8 3 x 4, 16 3 x 4 and for comparison the 
symmetrical lattice 16 4 . In general throughout these numerical simulations the masses of 
the different quarks were considerably heavier than the MILC97 data [50], which were taken 
between 7.5MeV < m q < 15MeV . The lighter quarks had masses in these newer simulations 
with the values between 40MeV < m q < 60MeV, while the heavier one was more in the range 
of lOOMeV < m q < 150Mel7, which are considerably higher than the light up and down 
quarks, but are not so far off for the strange quarks. The presented data was simulated using 
the p4 action for which comparisons were made to the more usual actions [72]. The resulting 
figures taken from these simulations we shall contrast the various different arrangements of 
flavors rif = 2, rif = 2 + 1 and rif = 3 given as functions of the temperature. In each of these 
flavor arrangements the quarks are at least an order of magnitude heavier than those in the 
previous figure 0] for the MILC97 data [50]. 

We now look into the properties of each thermodynamical quantity more specifi- 
cally. The Figure [7| shows the pressure p(T) for the three different flavors [70-72] in a way 
similar to that of the pure SU (3) gauge theory. Although the general shape of these curves are 
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Figure 7: The pressure for the three values of rif = 2, 2 + 1, 3 with massive dynamical quarks. 

quite similar in appearance, the values of the temperatures for the transitions are noticably 
lower in the theory with the dynamical fermions. Again the pressure was computed in the 
same manner as that of the pure lattice SU(N C ) using the integration method [8-10] as was 
discussed in the last part. In the case of the pressure p(T)[GeV/ fm 3 } the shape of all three 
curves is almost the same except for the starting point for the case of rif = 2. This is because of 
the difference in the value of T c . For the case of rif = 2 the value of T c = 175MeV is somewhat 
higher than rif = 2 + 1 and rif = 3 of around 155MeV. However, when one plots the pressure 
ratio p/T 4 , as was originally done [70-72], one notices a large difference in the different curves 
arising from the different number of degrees of freedom in each case, which changes the rate 
of approaching the Stefan-Boltzmann limit as well as the actual value of number itself. In our 
figure [7| the pressure in physical units represented on a logarithmetrical scale appear to start 
very near to the actual value of T c . However, in the form of the ratio p/T A the small values for 
the different flavors start well below T c with a very gradual increase until reaching the critical 
temperature. Thereabove the increase in each curve is at different rates due to the different 
number of degrees of freedom and the corresponding masses. In the actual original plot [72] 
(see figure 4.8) of p/T 4 against T/T c , one sees that the curve for rif = 2 lies well under that 
of rif = 2 + 1, which itself is significantly under that of the cases of rif = 3. Furthermore, the 
placing of these curves is quite different on the approach to the Stefan-Boltzmann limit. At 
the temperature of about 4T C the case of rif — 2 + 1 has only achieved around 75%, whereas 
both the cases rif = 2 and rif — 3 have arrived at about 80% of their respective limiting 
values for high temperatures. Whether this limiting value is actually attained in the high 
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temperature limit, still remains as an open question. Finally we remark that in the case of the 
pressure the errorbars are generally very small so that their presence is not essential to the figure. 

We next turn our attention to the energy density e(T) which presents quite a dif- 
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Figure 8: The energy density for the three values of nf = 2,2 + 1,3 with massive dynamical 
quarks. 

ferent problemmatic in the computation. The corresponding accuracy of the evaluations is 
considerably lower, which may be apparent by the presence of the errorbars in figure |H1 There 
we see that e(T)[GeV/ fm 3 } shows a quite different behavior as a function of the temperature 
T[MeV] from that of the pressure in the Figure [7| Even well below T c there are small values of 
s(T) which then grow much more rapidly at T c . In contrast the pressure really starts to take 
on sizable values only at temperatures very near to T c . Furthermore, one can see quite different 
rate of rise in e{T) depending upon the masses and the number of flavors. Nevertheless, its 
values start off at lower temperatures and rise more slowly than the corresponding values of 
p. Then just above T c the energy density s(T) rises very rapidly as a function of T[MeV]. We 
notice that the curve for n/ = 2 starts later because of the larger T c and rises more slowly 
just above T c . For temperatures above about 200MeV r e(T) for all the three flavors rise at 
approximately the same rate-practically on top of eachother. However, the nj = 2 the values 
always remain somewhat below the others as it was also the case for the pressure curves. 

From the knowledge of the pressure p(T) as well as the energy density e(T) we are able 
to calculate the equation of state in terms of e(T) — 3p(T). However, even as a difference 
between these basic physical quantities the equation of state still shows another type of 
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Figure 9: The equation of state is evaluated for the three values of nj = 2, 2 + 1, 3 including the 
massive dynamical quarks. 

behavior in relation to the flavors Uf and the quark masses m q . The values of the masses 
in the stated ranges have always held the nf = 2 values of these thermodynamical functions 
of the temperature to be smaller than nf = 2 + 1 and nf = 3. However, we see in figure HU 
that this statement holds only up to about 200MeV, above which temperature the nf = 2 
values remain within the errorbars clearly larger than the others. Furthermore, above about 
300MeV the case with n/ = 2 + 1 takes on larger values than rif = 3, but still remains smaller 
than nf = 2. Furthermore, in this range between 200MeV and 300MeV the change is very 
slow when compared with the region between 150MeV and 200MeV near the critical points. 
Thus we can clearly see here the types of contrasts between the parameters and the physical 
quantities. We remark also that for the equation of state the problem with the error bars at 
higher temperatures becomes very significant so that in all cases above 500MeV the curves 
give no real physical predictions. 

In this section we have numerically evaluated using the lattice gauge simulations [70-72]. 
We have shown the thermal properties in the three separate figures I7|8I and 01 the basic 
thermodynamical quantities- the pressure p(T), the energy density e(T) and the equation 
of state in the form e(T) — 3p(T) as functions of the temperature T. Furthermore, we have 
looked at these quantities in terms of the input parameters like the number of quark flavors 
and the quark masses, for which earlier in this part we have brought in the properties of chiral 
symmetry through the chiral condensate ratio. In the next section we shall briefly discuss 
some very recent data appearing within the last year, which we can compare with the above 
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shown results for different values of the lattice mass parameters. 



III. 5 Discussion of Recent Massive Quark Data 

Since the start of the actual writing of this work, there have appeared some newer 
data [73,74] for the three quark equation of state. In this section we shall present a discussion 
of some of these newer numerical simulations for 2+1 flavors with a brief analysis of the 
lattice results. In particular, we look at the simulations [74] plotted as a function of the 
ratio T/T c in their Figure 1, in which they include the data points for the following lattice 
quantities: the interaction measure, here written as I/T 4 , the pressure ratio p/T 4 and the 
energy density ratio e/T 4 . In their simulations they use a Symanzik improved gauge action 
and the Asqtad 0(a 2 ) improved staggered quark action for lattices with the temporal extents 
N T = 4 and N T = 6. They set their value of the heavy quark mass near to the physical 
strange quark mass m s . Then they choose the two degenerate light quark masses to have 
the values of 0.1m s and 0.2m s , for which they compute these quantities in the temperature 
range from 0.7T C to 2.15T C . For the computation of these thermodynamical ratios the integral 
method [8] has been used as discussed in previous sections. Furthermore, the estimated 
value [73] of the critical temperature T c is somewhat higher than for the Bielefeld 2+1 fla- 
vor data for this system with much lighter quarks, which is found to be about 1Q9 ±12 ±4MeV . 

These more recent results may be properly compared to the Bielefeld data [70-72] 
only for the simulations with N T = 4. In general we can compare the values for the interaction 
measure in the 2+1 flavor case [72] used in the above figure M to compute the equation of 
state to those values of I/T 4 plotted for N T = 4 in the Figure 1 of this newer work [74]. This 
comparison shows that the actual difference between the does not result in a large change 
in the interaction measure. Therefore, we could expect to have rather small changes in the 
properties of the equation of state over the temperature range in common to both cases. 
A simple pointwise comparison shows that for temperatures below and just above T c the 
values for the interaction measures are close to the same for both cases within the respective 
errorbars. Above 1.5T C the newer values of the interaction measure are actually somewhat 
larger. In the case of N T = 6 the numerical values for the interaction measure are considerably 
larger than the corresponding values for the Bielefeld data [70-72] just above T c and remain 
so throughout the higher temperatures. 

In summary for this new data we notice that the N T = 4 is generally closer to the 
numerical values used here. The larger values for the light quark mass provide the bigger 
numbers for the interaction measure especially around the critical temperature. We should 
also mention some values for the quantity e/T 4 approach the Stefan-Boltzmann limit for the 
N T — 6 simulations near to 1.25T C , after which the changes in value are quite small. This fact 
accounts for the larger growth of the interaction measure near to T c . This new data describes 
the QCD thermodynamics for three flavors of improved staggered quarks quite consistently 
with the data use here [70-72]. In the next part we will consider the results from our previous 
analyses of the equations of state in the figures E| and M in more detail in relation to the gluon 
and quark condensates at finite temperature. 
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IV. Gluon and Quark Condensates at Finite Temperature 



In this part we start our discussion of the thermal properties of the gluon condensates 
by using the approach that we described in the Introduction which includes both the pure 
SU(2) and SU(3) gauge theories [10,37] as well as the massive quarks [52]. We shall first 
discuss how the pure gluon condensate looks with no dynamical quarks present, for which the 
results for the equation of state in Part II can be directly used. The effects of the chiral phase 
transition with massive dynamical quarks are then discussed in various cases, for which we 
discuss the different relationships to the quark condensates. 



IV. 1 Pure Gluon Condensate 



The results for the gluon condensation in the cases of the pure gauge theories are shown 
in the figure El for the lattice data [9, 10]. In both these cases we have taken the zero 
temperature gluon condensate [36] to be {G 2 ) = O.O^fGeV 4 ] as the starting value for the 
gluon condensate at finite temperature, which we shall write for simplicity G 2 (T). We can 
see that this value remains very constant for temperatures up to nearly the deconfinement 
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Figure 10: We show the gluon condensates for the SU(2) and SU(3) invariant gauge theories, 
where both the ordinates are in [GeV 4 ]. 
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temperature [28,37], which are given as T d = 0.290 and 0.264 GeV for the color SU{2) and 
SU{3), respectively. The basic equation comes from the finite temperature trace anomaly 
which remains much the same as the one considered above used by Leutwyler [28] except that 
it is applied to the pure gauge theory not massless quarks. In the course of this section we 
shall reconsider the values in order to include the massive quarks. 

IV. 2 The Effect of Quark Condensates 

In the presence of massive quarks the trace of the energy-momentum tensor takes the 
altered form [44] from the trace anomaly as given by the equation © for which we recall 
that m q is the light (renormalized) quark mass and ip q , ip q represent the quark and antiquark 
fields respectively. In the last section we have discussed the temperature dependence of the 
chiral condensate in relatiion to chiral perturbation theory [28]. We recall that the second 
term in equation (ffilj) is quadratic in the temperature. This term corresponds to the process 
of pion absorption and emission given by the matrix element (it^ip^) . This term as well as 
the succeeding even powers in the expansion of {^pip) T allow the chiral condensate to melt 
quite gradually. Furthermore, we should note that in equation that all the terms subtract 
off the zero temperature condensate. In contrast to the multitude of processes in the chiral 
condensate melting, Leutwyler clearly points out [28] that such terms are absent in the pure 
gluonic processes since the field strength operators in the color combination given in equation 
(J3J) in the form G^G a is a chiral singlet so that the corresponding single pion matrix element 
vanishes [28]. When one calculates the corresponding trace of the energy momentum tensor, 
one finds for two massless flavors 

QHT) = ^-^-{ln^ - \\ + 0(T 10 ), (33) 
1 270 F£\ T 4 J v ; ' v ; 

where the logarithmic scale factor A p has the approximate value of 275MeV [28,45]. Thus 
in the framework of chiral perturbation theory the single and double loop contributions both 
vanish in the gluon condensate at finite temperatures. The first actual term present in the 
equation (J3~3~j) has the power T 8 arises from the three loop graph. The result of equation (J3llj) 
can be used in Leutwyler's equation (jHJ) to evaluate the finite temperature color averaged gluon 
condensate, which was discussed [28]. In the following we shall use these averages together in a 
single expression, both of which are evaluated from the numerical lattice gauge simulations with 
massive dynamical quarks. Furthermore, we shall include with these averages the renormal- 
ization group functions (3(g) and ~f(g,m), which appear in the trace of the energy momentum 
tensor even in the vacuum [46,47] from the renormalization process. 

Now we need to discuss further the changes in the computational procedure which arise 
from the presence of dynamical quarks with a finite mass. There have been recently a number 
of computations of the thermodynamical quantities in full QCD with two flavors of staggered 
quarks [48-50], and with four flavors [51,86]. We have already mentioned the problemmatic of 
the simulations for the dynamical quarks in the equation of state in the last section. Further- 
more, we mentioned there the problem of the quark condensate in relation to the restoration of 
chiral symmetry. Now we bring these two aspects of dynamical fermions together. Nevertheless, 
we first remark that these calculations are still not as accurate as those in pure gauge theory 
for several reasons. The first is the prohibitive cost of obtaining statistics similar to those 
obtained for pure lattice gauge theory. So the error on the interaction measure is considerably 
larger. The second reason, which is perhaps more serious, lies in the effect of the quark masses 
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currently simulated. They are still, in most of the cases in the present simulations relatively 
heavy, which unduely increases the contribution of the average quark condensate part in the 
interaction measure. In fact, it is known that the vacuum expectation values for very heavy 
quarks is proportional to the gluon condensate Gq, written in simplified notation, which in the 
first approximation is given by [36] 

($M = j^-Gl- (34) 

Furthermore, there is an additional difficulty in setting properly the temperature scale even 
to the extent of rather large changes in the critical temperature have been reported in the 
literature depending upon the method of extraction. For two flavors of quarks the values of 
T c lie between 0.i40Ge^ [50] and about O.lTOGeV^ [63] which is considered presently a good 
estimate of the physical value for the critical temperature. 

We are now able to write down an equation for the temperature dependence of the 
thermally averaged trace of the energy momentum tensor including the effects of the light 
quarks with a mass m from A m (T) so that 

6^(T) = A m (T)-T 4 . (35) 

The thermally averaged gluon condensate is computed including the light quarks in the trace 
anomaly using the equation (JBJ and the interaction measure in 0p(T) to get 

G 2 (T) + J2 q m SM T = G l + E ? m «^>o - ^( T )> ( 36 ) 

It is possible to see from this equation that at very low temperatures the additional contribution 
to the temperature dependence of the gluon condensate from the quark condensate is rather 
insignificant and disappears completely at zero temperature. However, in the range where the 
chiral symmetry is being restored there is an additional effect from the term (ijj q ijj q )T, which 
lowers G 2 (T). Well above T c after the chiral symmetry has been mostly restored the only 
remaining effect of the quark condensate is that of m q {^) q ij) q )Q. It is then known [36] how this 
term shifts the gluon condensate of the vacuum. Thus we may well expect [37] that for the 
light quarks the temperature dependence can only be important below and near to T c . In the 
case of the chiral limit m q — > the equation (|37|) takes the form of Leutwyler's equation © 
as, of course, it should because Leutwyler used two massless quarks [28]. For the smaller values 
of the simulated quark masses in lattice units of 0.01 to 0.02 ('ipqip^T has mostly disappeared 
in the range where G 2 (T) differs from (G 2 )q. 

IV. 3 Gluon and Quark Condensates with two light Flavors 

In the way we mentioned above for equation (j3J) with the change of the pure gluon 
condensate to the sum of the quark contributions (G 2 ) + Y^ q m q('4 , q l l ) q) ■ First we discuss the 
finite temperature gluon condensate in the "chiral limit" m q — > as shown in the figure ITT1 
with the two m q = curves. From the MILC light quark data [50] we are able to subtract the 
quark condensates at finite temperatures multiplied by the two quark masses m q = 7.5MeV 
and m q = 15MeV for the two sets of lattice data. The effects of the temperature as well as the 
restoration of the chiral symmetry are included with the decreasing value of m q < ipip > with 
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increasing temperature are gotten from the original lattice data [50]. We use for the vacuum 
value of < ipip >o= —(259 ± 27MeV) 3 from a recent estimate [75] taken from newer MILC 
dynamical quark data [76]. Now we can reformulate the above equation (|36|) by replacing the 
sum over the flavors simply by the number two. The the finite temperature gluon condensate 
is given by 

G 2 {T) = G 2 + 2mSMo (l - ®M) - ^(T). (37) 

The ratio of the finite temperature quark condensate to its vacuum value, (V^V^r/ (VyA?)o> 
we write simply as (ipqipq), which is just the chiral order parameter. The values can be seen 
in the figure 4 of the MILC data [50]. We show these results in the figure ITT1 for the two 
different stated masses including the errorbars arising from both the interaction measure and 
the quark condensate data. We notice that with increasing temperatures that the finite quark 
mass lowers the gluon condensate. Furthermore, it can be seen in figure ITT1 that between the 
temperatures of about 140MeV and 160MeV, which is just below T c , both with and without 
the quark condensates lie very mush together- almost within mutual errorbars. Below and 
above these values the heavier quark values appear to digress although there are simply too 
few points to allow for a firm conclusion. This fact clearly shows an additional effect from the 
restoration of the chiral symmetry on the thermal gluon condensate in the presence of massive 
quarks. At the highest temperature around 240MeV the chiral symmetry for the smaller 
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Figure 11: The finite temperature gluon condensate is shown for two light quarks given by the 
two different lattice masses. We compare these with the corresponding data in the chiral limit 
shown with the triangles. 
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mass value is about 95% restored. From these results we can understand how the effect of 
the quark masses together with the restoration of the chiral symmetry changes the size of 
the gluon condensate. Whereupon we may conclude this discussion of the gluon condensate 
with two light quark from the MILC97 data [50] by saying that the decondensation of gluons 
begins considerably below the critical temperature and then slows down for a temperature 
range above it. Finally we remark that we have taken a compromise value for the vacuum 
expectation value G$ = 1.95[GeV/ fm 3 ]. This choice for the simulations with dynamical quarks 
is somewhat larger than the earlier value [36] used for the pure gauge theory in the figure ITU1 
However, it is well within the range of many of the later estimates (for the literature see [15,16]). 

IV. 4 Gluon and Quark Condensates with more massive Flavors 

After this evaluation of the gluon condensate for the very light quarks we present the 
results for two and three moderately light quarks using the Bielefeld data [70-72]. Thereafter 
we look at the case with two lighter quarks together with one heavier quark. We shall in 
general compare the effects of the pure gluon condensate with the maximally estimated effects 
arising from the quark mass terms, from which we can derive the changes in the thermal 
properties of the gluon condensate which are due to the dynamical quarks. 

We see in Fig. H2] the effects of decondensation of colored quark-gluon fields containing two 
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Figure 12: The pure gluon condensate < G 2 >t in units of GeV/ fm 3 for two moderately light 
quarks is indicated by the vertical/horizontal lines. The crosses show the estimated maximal 
effects of the (thermal) quark masses. 

fairly light dynamical quarks. Here we can clearly notice that above 200 MeV the difference 
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between the pure gluon condensate with the errorbars and the mass terms appearing with the 
chiral restoration given by the crosses. The darkened region between the horizontal line and 
the cross x is the immediate effect of the thermal mass term caused by the chiral restoration. 
We can see here that the thermal effects of the mass terms coming from the quark condensate 
do partially slow down the decondensation of the gluons. This effect we can compare with one 
of our earlier analyses [64,65] of the data for nj = 2 using the even earlier MILC data [48-50]. 
For the MILC data the quarks were much lighter- almost an order of magnitude less- from 
less than 5MeV to around 15MeV. We can now compare the numerical data in the last 
section for the equation of state shown for different flavor numbers in the figure 03 The gluon 
condensate we shall denote here as (G 2 )t- In the following figures we shall leave out the data 
for the temperatures above 500MeV because of the very large errorbars. Unfortunately, in 
the following analyses we have lacked the comparable data for the chiral averages for the 
quark condensate (■0 ? ^)t which may be associated with the temperature as was the case 
for the MILC data [50]. Thus the shaded region for the curve above T c gives the estimated 
maximal range of the effect of the chiral condensate, which in some cases may reach beyond 
the errorbars. 
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Figure 13: The gluon condensate< G 2 >t for three middle light quarks in physical units 
[GeV/fm 3 ]. 

Next we show in the figure the effects of three moderately light quark flavors on the 
gluon condensate (G 2 )t- A quick examination of this plot shows that the change of (G 2 )t is 
not so rapid over the same temperature interval as it was the case for the two quark flavors 
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with similar values for the masses. Again we also show here the maximal estimated effect 
coming out of the masses which when it is taken into account raise the points in the chiral 
limit. This same effect we have already just seen in the MILC data in figure El The small 
cross x marks the estimated highest value for (G 2 )t due to the chiral restoration. 

We show in the next figure El the computed results due to the effects of two lighter 
quark flavors together with one much heavier quark. A quick examination of this drawing 
shows that the large masses of the dynamical quarks flattens the descent of the decondensation 
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Figure 14: The gluon condensate< G 2 >t for two moderately light quarks and one much heavier 
one. 

curve even further in comparison to the other two cases. Clearly in this case the mass effects 
are much larger, which is primarily due to the much more massive "strange" quark. We 
see that at higher temperatures the crosses x are generally above the errorbars for the pure 
gluon condensate. It is obvious in this case with a heavier quark that we are very far from 
the situation of the chiral limit. Thus the restoration of the chiral symmetry has much less 
effect on the gluon condensate in the temperature region above the critical temperature T c 
and below the deconfinement temperature of the pure SU(3) gauge theory. This remark is 
quite consistent with what we had noted in the last part in the figure El where we compared 
the effects of the quark masses on < ifjip >. We saw there that the larger mass value caused 
the curves for < ipip > to move towards higher values in the coupling (3. Furthermore, the 
descent of the curves at higher masses is slower, which means that the resoration is less rapid 
at T c for the more massive dynamical quarks. 
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IV. 5 Comparisons of Properties of Gluon Condensates 

In this segment we discuss some previous work [65] on gluon condensates. In particu- 
lar, we shall take some special extreme examples out of the cases for which we have numerical 
and analytical evaluations. In the previous part of this report we have presented some 
numerical results for the gluon condensates- first for the pure lattice gauge theory, then for the 
light dynamical quarks with two flavors and finally for the somewhat heavier dynamical quark 
flavor combination. Therein we had also made some comparitive remarks between the various 
cases. However, in the preceeding part we have not tried to carefully compare the results 
for each system. Some years ago such a comparison had been discussed between the gluon 
condensates arising from the pure SU (3) lattice gauge simulations [10, 37] with those coming 
from the light dynamical quark data from the MILC collaboration [50]. As an additional 
point of comparison [65] a gluon condensate arising from a pure ideal gas of gluons was also 
introduced with the same ground state value as the QCD vacuum gluon condensate Gq, for 
which the 0.012GeV A in often taken [36]. For the ideal gluon condensate the condensation 
temperature T was chosen to be T c , the deconfinement or critical temperature of the pure 
lattice gauge theory. We show the results of these comparisons in the figure 

We now briefly discuss the condensation of an ideal gas of gluons, which we simply write as 
G 2 (T) for a finite temperature T. The vacuum expectation value of the pure gluonic system we 
denote by Gq, which has the above known value at zero temperature [36]. From dimensional 
considerations of the structure of the gluon condensate we can write down the ideal gluon 
condensate for T < T 



where T is the temperature at which the condensation process ceases. For T > T Q then 
G 2 (T) = 0. This means that for the gluon condensate of an ideal gas system in equilibrium 
with its groundstate contribution must vanish above the condensation temperature. This form 
is just a generalization of the relativistic Bose-Einstein condensation to a four dimensional 
Euclidean space with T acting as the critical temperature. Therefore, this equation gives 
a simple relation between the gluon condensate at a finite temperature T and that at zero 
temperature. Thus we relate Gq to the QCD vacuum condensate and To to the deconfinement 
temperature, which are the corresponding quantities in the previous parts. Although this 
form for G 2 (T) as an ideal gluon gas seems at the present time very highly oversimplified, 
many of the earlier analyses in QCD at finite temperatures with gas models had assumed free 
gluons. This includes some earlier collaborative work of the author [77] where this type of 
thermal behavior was assumed as a first approximation to the decondensation of gluons. The 
extent to which this assumption is valid we are now able to evaluate from the pure SU(3) 
data [10]. Furthermore, we will make no attempts here at the inclusion of the hadrons in the 
thermodynamics of the gluon condensate. 

At the beginning of the previous part we have shown the results [37] for the gluon 
condensates of the SU(2) and the SU(3) invariant lattice gauge theories over the whole 
temperature range up to OAGeV. In both cases the uppermost temperature is well above 
the deconfinement temperature, for which the values of G 2 (T) are all less than zero as is 
shown in the figure In the following paragraph we used the argument of Leutwyler [28] in 
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Figure 15: The lines show the gluon condensates for SU(3) (solid) and the ideal gluon gas 
(broken) in comparison with that of the light dynamical quarks denoted by the open circles 
and the heavier ones with filled squares. The error bars are included when significant. The 
physical units are for G 2 {T) [GeV 4 ] and T [Gel/]. 



order to explain why the pure gluon condensate remains constant while retaining its vacuum 
value Gq almost until the temperature 9 reaches T c . This effect clearly comes from the strong 
QCD attractive forces between the gluons. Just below T c the drop in G 2 {T) signifies the 
start of the decondensation process. Not only do the gluons coming from the interactions 
in the vacuum condensate decondense, but the further gluons which are created at these 
high temperatures also take part in the decondensation process. We now can understand the 
difference in the behavior for the condensate of the ideal gluon gas which is shown in the 
figure ED This gas of free gluons shows a continual decrease even at the low temperatures well 
under T c . At and above T c the ideal gluon gas condensate vanishes so that G 2 (T) = for all 
T >Tq. However, the pure SU(N C ) gluon condensate continues into negative values of G 2 (T) 
as it was shown before in the Figure Thus the difference between the pure lattice gauge 
thermodynamics and the ideal gluon gas is apparent over the whole range of finite temperatures. 

9 H. Leutwyler, "Quark and Gluon Condensates", plenary lecture at QUARK MATTER '96, Heidelberg, 
Germany on 24. May 1996. Some further details on this approach were given by Heinrich Leutwyler on this 
subject in a private discussion with Graham Boyd and the author [37] at Heidelberg. Here we take T c as the 
critical temperature for the pure gluon gas. 
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The comparison of the gluon condensate from the two light dynamical quarks brings 
new issues into these considerations. As we have discussed in the preceeding two parts 
of our work, the presence of the chiral restoration transition changes the thermodynamics 
completely. In the equation of state e — 3p for the case of the two light quarks we can 
see in figure 0] how around 140MeV both the curves rise very rapidly. However, at about 
200MeV the lighter quark mass slows up and even appears to drop a little. Not only is 
this change at much lower temperatures than the pure gauge results in the figure ^ which 
show a steady rise beginning at a much higher temperature. The presence of the errorbars 
in the figure |U show that the apparent decline of the last point at about 240MeV is well 
within the errorbars of the preceeding point at about 196MeV. The MILC data is then 
again plotted for the gluon condensate G 2 (T) in the figure ITT1 for which the units are taken 
as [GeV/ fm 3 ] in contrast to the above figure El In figure ^2 we have connected the points 
together in order to provide a better view of the changes. The two curves with the errorbars 
have included the effects of the masses with the quark condensates, which adds some further 
errors. Nevertheless, we see that the gluon condensate then continues to fall slightly. The 
unconnected points in the above figure El correspond to the unconnected points given in the 
chiral limit. This fact has been recently pointed out [24,25] in relation to the "soft" glue, 
which disappears much earlier than the hard glue or "epoxy". It can also be related to the 
problem of mass and the presence of mesonic bound states [29] above T c , which here denotes 
the actual critical temperature for the two flavor light quark thermodynamical system. We 
shall return to this topic in relation to chiral resoration models in a later section of the next part. 

Finally we are now in a position to discuss the interrelationship between the process 
of the gluon decondensation and the restoration of chiral symmetry. In the previous three 
figures El El and E] we can see the varying extents that the dynamical quark properties 
affect the amount the two different broken symmetries are altered with the change of the 
temperature. As it has been clearly stated by Leutwyler [28], the quantity < ipip > acts 
like an order parameter for the breaking of the chiral symmetry. Actually it is the term 
1 — < ipip > in the equation (J37j) which acts as the restoration quantity for the effects 
of the quark condensate on the gluon decondensation. In the figure ^2 we can see the 
direct effects of this term which clearly lowers the gluon condensation curve with increasing 
temperatures due to the restoration of the chiral symmertry. Naturally the bigger effect 
of the decondensation of the gluons arises in the interaction measure A m (T) in equation 
(|26|) from the mass renormalization of the quarks. We see that the antiquark-quark average 
over the lattice or site average clearly lowers the plaquette average from the links in the 
interaction measure. This we had already seen in figure 0] for the equation of state for 
the light dynamical quarks. The equation of state for the dynamical quarks still provides 
the most important contribution to the decondensation of the gluons. However, the extent 
to which this last statement is true varies with the number of flavors and the masses 
of the quarks. Thus we can clearly see the main points of contrast between the deconden- 
sation of the gluons in the pure gauge theory and that of the theory with the dynamical quarks. 



36 



V. Discussion of Physical Results 



In this part of our work we shall describe how we are able to calculate some physical 
quantities from the structure of the equation of state which relate to such phenomena as the 
possible phase transitions between the hadronic states and the much sought after quark-gluon 
plasma. First we discuss some models relating to the relativistic high energy collisions at high 
temperatures and large energy denstiies [15]. Thereafter we shall look at the low energy and 
low temperature limits where the effective Lagrangians are very useful. We discuss a model for 
scalar meson dominance with the coupling to the trace of the energy-momentum tensor [18]. 
In the last section of this part we look at some other aspects of the equation of state which 
arise at temperatures well below the critical temperature. In our last comments at the end we 
mention how these considerations may relate to some newer lattice simulations not included 
in the presented data. 



V.l Phenomenological Models of Quark and Gluon Properties 

One of the most successful models for the highly relativistic nucleus-nucleus collisions 
was proposed almost a quarter of a century ago by Bjorken [91]. He suggested a senario 
where a quark-gluon plasma should survive on a time scale of 5/m/c, where c is the speed 
of light which we have usually taken to be unity. During this short time the collision region 
maintains a longitudinal flow which should be described according to the laws of relativistic 
hydrodynamics [92]. 

The dynamical structure of the Bjorken model can be reduced to that of a 1+1 di- 
mensional Lorentz boost invariant system in the time and space coordinates t and z. Bjorken, 
however, then chose to represent the flow in terms of the "light cone" variables, for which he 
took the propertime r = y/t 2 — z 2 and the rapidity y = \ In These variables are useful in 
maintaining the boost invariance. In addition to the usual energy momentum conservation 
Bjorken further postulated [4,15] a conservation of the entropy density s = (3{e + p), where 
(5 is the inverse temperature. Out of these assumptions one gets for the fourvector entropy 
density s M = su^ 1 the conservation law 

= 0, (39) 

where VtF" is the fluid fourvelocity with u^u^ = 1. After setting the initial propertime tq with 
the entropy density s(r ), one finds simply the solution [91] 

44 = X (40) 

In terms of the variables r and y the entropy density stays on the lightcone. This result 
may be interpreted as the entropy content per unit rapidity is a constant of the motion [91]. 
Furthermore, this result may be related to the Bjorken equation of state, which may be written 

as 

e=p(3 + A 1 ). (41) 
Comparably the entropy density for the Bjorken model may be written as 

s = pp(4 + A 1 ). (42) 
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The new quantity Ai is denned in the Bjorken model as follows: 



T ( dn . 

A i = ~ 3^ ■ ( 4 3) 



n \dT 

In the rest of this very original work Bjorken [91] showed how these equations can be used 
for the effective number of plasma degrees of freedom n(T) which he had originally proposed 
to have a sharp rise for the temperatures above the transition temperature around 200MeV. 
An important application of an extension of the Bjorken equation of state (}4*2*|) has been found 
in the calculation of the radiative energy loss of high energy partons traversing an expanding 
QCD plasma [93]. 



For us in this work what is really significant here is that the Bjorken equation of 
state leads to corrections beyond the ideal gas as well as the simple bag model pictures. 
Bjorken [91] also pointed out that the general features of the hydrodynamical expansion 
follow from the positivity condition on the trace of the energy momentum tensor. In fact, 
this statement would imply that the positive definiteness condition on the trace of the energy 
momentum tensor is such that the trace itself would only be zero in the absence of all 
interactions. 



More recently the hydrodynamical models [15] have been brought in relation to the 
Bjorken model. In this context the well known ideas of relativistic hydrodynamics [92] have 
been extended for use in the explanation of the data from heavy ion collisions undergoing 
a hydrodynamical expansion [94]. At the present time there is a considerable amount of 
phenomonological work in this direction. 

V.2 Theoretical Models for the Quark and Gluon Properties 

In this section we shall look into some of the theoretically important statistical models 
for strong interaction physics. Some of these relate to the recent RHIC data [15]. First we 
mention the Statistical Bootstrap Model (SBM) which was earlier regarded for providing the 
limiting temperature for hadronic physics [87]. In the following paragraphs we shall discuss 
here some further models involving the quarks and gluons in large statistical systems. 



Now we discuss the properties of resonances at high temperatures and densities which are 
typical of high energy particle collisions. The particle creation process at such high energies 
is the starting point for the SBM first introduced in the midsixties by Rolf Hagedorn [87]. 
Earlier pure statistical models based on ideal gas theory with the proper statistics related 
the temperature directly to the center of mass energy. These models had in common an 
exponential energy distribution relating to the probabilities of each of the hadronic states. 
However, Hagedorn had noted that in high energy proton-proton collisions that the masses m 
are distributed in what appeared to be asymptotically an exponentially rising mass spectrum 
p(m). This spectrum included the entire known resonance structure of hadronic particle 
physics. Furthermore, he then postulated that the mass spectrum approached asymptotically 
the energy level density cr{E) in the high energy limit. After having calculated c(E) from the 
partitiion function with p(m) for a gas of particles and resonances, he invoked the statistical 
bootstrap condition relating p(m) asymptotically to cr{E). However, this condition only held 
true up to a particular temperature T = 158 ± 3[MeK], when certain statistical assumptions 
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were made on the pararmeters of p(m). Nevertheless, this temperature T had found its 
special place in the particle phenomenology of that time through its place as the limiting 
temperature. A more contemporary understanding of this resonance formation is found from 
an examination of this Hagedorn temperature and the partition thermodynamics [88]. The 
hadronic phase with all the included resonances can only exist below T . This statistical theory 
was clearly modeled on the particle states which were known at that time. The sense of such 
a model with a limiting temperature has changed drastically with the inclusion of quarks as 
the basic constituents of the hadrons. More recently some further implications in relation to 
the quark-gluon transition have been discussed [14]. The SBM serves as a starting point of a 
statistical treatment of the hadrons with their constituent quarks and interacting gluons. 

In the early nineteen seventies not long after the initial development [16, 20, 56, 57] 
of QCD as the gauge theory for the strong interactions between the quarks, there was a great 
deal interest in finding analytical solutions for particular cases. The discovery of the instanton 
solution [79] in 1975 of the SU(2) Yang-Mills field and its further development [15, 17] were 
important starting point for an approximate analytical approach to nonperturbative QCD. 
After a brief initial excitement giving hope of a discovered particle in the years directly 
following these discoveries, there began an era of long and systemmatic investigation of these 
rather intriguing mathematical objects [80,81] in the name of the instanton gases and liquids 
with various forms [15]. More recently the importance of the bound states and the instanton 
molecules [83] has become clearer. Furthermore, it was noted that the fermionic (quark) zero 
modes can be related to the axial anomaly [21]. This leads to the properties of the chiral 
symmetry breaking and its restoration [25] as we had mentioned previously. A more lengthy 
discussion of these ideas can be found [15,16] which give details on the explicit structure. The 
relation to some earlier lattice computations is also quite important (see [5] in section 17.6). 
Furthermore, the thermodyanmical properties of the instantons relating to the topological 
charge Q\ can be evaluated on the lattice using the thermal susceptibility in equation (J3T|) 
as indicated above in the equation (J3~2~|) . There are also related lattice simulations on the 
flux tubes (see [5] section 17.7) for further discussions on the newer evaluations. Further 
work involving the instanton has been presented by Ilgenfritz and Shuryak [82] in a simple 
model at finite temperature. The role of the instantons has been generally discussed by 
Shuryak [15] in relation to the transition from hadrons to the quark gluon plasma. As we 
mentioned in the last part, there was an attempt [77] to formulate a mean field model of an in- 
stanton (caloron) gas at finite temperature in contact with the condensate of an ideal gluon gas. 

The early finite temperature simulations [2, 5] for lattice gauge theory provided in- 
formation of the thermodynamical functions around the critical temperature T c . The 
discussion of some previous results for the gluon and quark condensates were derived from 
some of the earlier simulations. In this relationship we must mention that the values of the 
vacuum condensates were calculated out of sum rules [36,96]. The early use of the finite 
temperature lattice results [22, 27] has had a significant part in thedevelopment in such 
computations. These authors had evaluated the gluon condensate in the separate parts arising 
from the chromoelectric and chromomagnetic fields with the thermally averaged condensates 
as seen in reference [27] in their Figure 4. Unfortunately in the earlier lattice data the 
temperature scale was very hard to determine. In 1994 this changed significantly through a 
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better evaluation of the lattice beta function [9] for the SU(2) gauge theory 10 . Even in these 
earlier times the use of these ideas was apparent in the study of the hadron to quark-gluon 
transition [14]. 

Finally we mention some other work involving QCD at finite temperature and den- 
sity. First we remark on the models involving quasiparticles [89, 90] which take into 
consideration the lattice properties of the high temperature states. These models generally 
reproduce the lattice data above about 2T C . Some recent lattice work [69] using a QCD model 
with two flavors of adjoint quarks provides a very unique picture of both the deconfinement 
and chiral phase transitions. It shows both a chiral second order transition and a first order 
deconfinement phase transition. However, it can be seen in their Figure 12 showing the chiral 
condensate in terms of the coupling j3 and the lattice quark mass m q a that in this model the 
deconfinement transition appears at a lower temperature than that of the chiral restoration. In 
fact they estimate [69] that T c /Td pa 7.8(2), which means that the critical temperature remains 
very far above the deconfinement temperature. 

In this section we have related the results of numerical simulations which are usu- 
ally represented as ratios of the quantities to powers of the temperature to the actual physical 
quantities like the pressure, energy density, entropy density. We can arrive at a form of the 
equation of state which can be used to find other physical quantities. We look again at the 
first pure gauge theory simulations at finite temperature, which provided for us a way to arrive 
at the thermodynamical quantities in the gluon condensate. 



V.3 Scaling Properties in Matter 



There are a number prominent examples where the parameters in the effective actions 
for physical processes need be determined by some further rules. Important among these 
is the Brown-Rho Scaling properties for finite densities p, which we have mentioned in the 
introduction [23,24]. It is commonly given by a scaling function for the different hadrons 
(for a general description see reference [24]) in the following form: 

*<p) » K <id « !I&M „ MliEl. (44) 

U m * fv M 

The expressions f*(p) and f n are the pion decay constants in the medium with the density p 
and that of the vacuum respectively. In this same notation (with the * above) the masses of the 
scalar a and vector V particles are compared with the corresponding values in their different 
media. In the last ratio there appears the nucleon mass M. Nevertheless, it is important in all 
cases to stress that the Brown-Rho Scaling is a mean field relation which emerges at the tree 
level from the effective Lagrangian for the chiral fields [23,24]. This scaling rule may be further 
approximately generalized to the quark condensates are related to the hadron masses. 

- - ¥t- (45) 
m (qq) 

10 The author the made an analysis of the trace anomaly relating to the lattice data from the interaction 
measure to the gluon condensate, under the title "The Trace of the Energy Momentum Tensor and the Lattice 
Interaction Measure at Finite Temperature", Bielefeld Preprint, BI-TP 94/41. This first calculation was ex- 
tended in later works [37,52,64,65] using the same approach in which data for SU(3) pure gauge theory as well 
as for dynamical quarks were included. 
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This case is referred to as the "Nambu scaling" which usually includes the temperature 
dependence [24,27]. This property is a generic feature of the linear sigma model. It is then 
supposed to be upheld as the temperature rises towards the chiral restoration temperature. 
It was noted in a recent work how the matter induced modifications of certain resonances 
fitting the masses of various particles [78] causes a shift as well as a change of shape of various 
particles and resonances. Further work on the nature of the chiral restoration [25], used an 
earlier work [65] of the author to discuss how the mixed phase [29] appeared after the chiral 
restoration has started, but before the deconfinement temperature in the pure gauge theory. 
We have already seen in the last part that the Figure E3 can be interpreted in terms of the 
soft glue, which disappears rather quickly with the chiral restoration, while the hard glue 
(epoxy) remains well above T c . There these properties were associated with the slowing up of 
the decondensation processes of the light dynamical quarks. The heavier quarks gave a much 
flatter decondensation curve for which the two different properties are harder to recognize. 
Furthermore, it is expected that the correlations above T c are, indeed, very significant and 
have an important relationship to the recent RHIC experements [24, 25] . 

V.4 Scalar meson dominance 

In this section we look into the model of Freund and Nambu [18] for the dominance of 
the scalar meson 11 . This property is deeply rooted in the empirical observations of the meson 
and baryon mass splittings. Already in 1968 they had proposed that the classical meson field 
(f(x) be coupled to the trace of the energy momentum tensor 9£. A more detailed discussion 
of the relation of these (pseudo) scalar fields to the energy momentum tensor one finds in 
section three of Bogoliubov and Shirkov [84]. In this model the nonhomogeneous Klein-Gordon 
equation could be written in the form 

(□ - m 2 Mx) = gB% (46) 

whereby m is the scalar meson mass and g its coupling to the trace. Here it is important 
for us to note that this coupling g differs from the usual QCD coupling in its place in the 
model, which was, of course, not known at that time. In the following paragraphs we briefly 
sketch the basic arguments for this model which we shall later relate phenomenologially to the 
process of gluon decondensation. 

From the usual properties of field theory one can derive the trace of the energy mo- 
mentum tensor [84] ^ from the Lagrangian C((p, d^cp). Out of the above form of the 
Klein-Gordon equation (J46j) for the scalar field (p(x) Freund and Nambu derived [18] their 
effective Lagrangian by taking into account the appropriate boundry conditions. It took the 
following form: 

£{g,cp) = h^R(x)-% V - \m\\ (47) 

where R(x) is defined to be (1 + 2g(p(x)). After a careful consideration of the form of the 
Lagrangian they [18] introduced a new field if)(x), which is defined by 

$(x) = g- 1 R{x) 1/2 . (48) 

11 The actual scalar meson states are discussed in the "Review of Particle Physics" [95] in a "Note on Scalar 
Mesons" pp. 506 -510 and 522 -526. The known low lying scalar meson states go under the names /o(600) (or 
a) and /o(980) for the isoscalar as well as ao(980) for the isovector. These states are all known to decay into 
various pairs of pseudoscalar mesons and, of course, secondarily into pairs of photons. 
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This field transformation yielded a new Lagrangian £'(g,ifj) in terms of the new field in the 
form 

1 1 1 m 2 

= + - ^W- gj. (49) 

Furthermore, it must be noted here that the sign is wrong for the usual "mass term" in the 
Lagrangian density. This situation provides the conditions [85] of a degenerate vacuum which 
gives rise to Goldstone bosons upon quantization. It can be easily seen after one further change 
to ip± as ip±l/g. Then we find this change of the fields results in 

£'(<?, V±) = \{d^ ± f - \m\l ± \m 2 g^ ± - \m 2 g 2 ^ (50) 

which is equal [18] to the original Lagrangian C(g,(p). As we would expect, the dynamical 
properties of both of these fields ip± correspond to the same mass m. Either one of the 
Lagrangians C(g,i[)) or C\g 1 ip±) yields after quantization the same physical properties in the 
corresponding S Matrix as the original one C(g,ip) derived [18] in equation (f4Tj) . Thereby 
it is clear that there is a twofold degeneracy in the vacuum of the quantized theory. This 
model even in its unquantized form was recognized by Bruno Zumino to break the scale 
invariance [85] as well as allow for the Goldstone bosons. Furthermore, this situation leads to 
the fact that the coupling of the trace g9£ in the nonhomogeneous term of the Klein-Gordon 
equation (|4fij) relates directly to the square of the mass times the strength of the transformed 
scalar field. Thus the Lagrangian of Freund and Nambu provides a simple model for the scalar 
field dominance of the trace of the energy momentum tensor 

The quantized theory can be achieved from the classical one in the usual way from 
the Poisson brackets [84]. The momentum tt(x) canonically conjugate to <p{x)) from the 
original Lagrangian C(g,<p) is gotten as the Poisson bracket ^{R(x)~ x , dof}- Because of this 
form of the Poisson brackets the ordinary equal-time commutation relations do not specify the 
canonical commutator but only the mixed commutator in the form 

\{Rtf)-\[<p(x),&w{2>)]}^ =a , Q = i5 3 (x-x'). (51) 

From this form and the specification that [<p(x), <p(x')] Xo=x ' and [da<p(x), d' Q (p(x')] Xo=x > o both 
identically vanish allows the calculation of the equal time field momenta [vr(x), ff(x')] XQ=x ' o . 
With the definition that at the same point 

[^(x)- 1 ,^)- 1 ]^) = 0. (52) 

it is possible to evaluate the Schwinger commutator [34,35] for 9oo(x). More important to 
us in the present work is the fact that Freund and Nambu have evaluated the corresponding 
commutator for the trace of the energy momentum tensor 0£(x) as 

[W^OOU^o = - 4 *(z) + 9o l (x')]d l 5 3 (x-x'). (53) 

This feature allows that the trace of the energy momentum tensor be treated as a quantum 
field. These authors [18] have suggested that this commutation relationship be considered as 
a general feature of the theory of hadrons. Thus the Lagrangian C(g, (p) in their formulation 
provides a rather simple model for scalar field dominance from the trace of the energy 
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momentum tensor. These results can be understood from the transformation properties of the 
action integral under the scale transformations [85]. 

Next we work out the special case of the scalar meson field (p(x) with the mass m 
for the ground state structure in the simple bag model. For simplicity we assume a simple 
cubic region with the sides of length l[/m] for the discrete momentum representation [84]. 
The solution for the zero mode of the nonhomogeneous Klein-Gordon equation (j4T)j) for the 
scalar field tpo{x) using the isotropic bag model value [14] of the energy momentum tensor 
6^ = 4£>. In this case we find a simple time dependent solution (p(xo) with the zero spatial 
modes for the scalar meson ground state solution 

^o(^o) = (cos(mx ) - 1) , (54) 
m 

where we have taken <^o(0) = and y?o(0) = 0. This special case is so easily solvable since there 
is no spatial dependence coming into the field equation from the trace of the energy momentum 
tensor. Thereupon the value of the amplitude only involves the nonhomogeneous term divided 
by the mass squared. Then this special case yields just the simple harmonic oscillator type of 
solution. Now we see for the chiral limit for the solution in equation (p>4*)) that the solution goes 
to the finite value 

limm — > 0, <Po(xq) — > 2gBx 2 . (55) 

Thereby the small mass limit for the zero momentum state has a finite value which is positive 
definite from xq 2 . Since this zero momentum state solution has no spatial dependence in it, we 
may replace xq by the propertime r at x = 0. Then we find at the center for limm — > that 

Mr) = 2gBr 2 . (56) 

Thus we would expect the scalar field for the hadron to grow quadratically in the proper 
time, for which only on the light cone does it totally disappear. Furthermore, after a little 
further investigation, we find that all the higher mode solutions in the discrete momentum 
representation would involve free spatial wave solutions set to zero at the boundries of the 
box. However, these solutions involve an infinite sum in each spatial dimension. 

As a concluding point to this discussion we comment on the interesting possibilities 
of this model from the view of the equation of state. Next we shall consider the trace of 
the total energy momentum tensor Tjf(T) as a function of the temperature. In the simple 
solution found above for T = we saw that the trace appeared only in the solution form the 
bag constant 4£>. Then we would expect that even at finite temperatures this same type of 
dependence on the equation of state would directly enter into the solution. For the sake of 
simplicity we shall now only consider the effects of the decondensation for the pure SU(3) 
gluons as shown in the Figure ITU1 We recall that the decondensation of gluons starts around 
the deconfinement temperature of about 270MeV, where the equation of state undergoes 
a rapid growth as shown in the first figure for the temperature dependence. If we assume 
for simplicity that the scalar meson state in the simplest case is only changed by the added 
temperature dependence of the trace, then we would expect a rapid growth in the scalar field 
<Po(x ) as a function of the temperature. Thus our boundstate solution for the nonhomogenous 
Klein-Gordon Equation (J4lj|) given in ([54|) . shows rapid changes in the classical scalar meson 
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field just above the decondensation temperature 12 . Furthermore, if we take the lowest scalar 
meson, the a, with a mean mass m a = 600MeV and a full width at least AOOMeV, this added 
energy from the decondensing gluons could perhaps be enough to bring about a single decay of 
a virtual scalar meson state into two pions. In no sense are we advocating that this very rough 
picture should provide a realistic explanation for the effects of gluon condensation in QCD. 
However, it is generally consistent with the decay of the scalar meson in its very simplest case. 

V.5 Entropy for the Hadronic Ground State 

In this section we look at a different situation involving the equation of state for the 
quarks making up the ground state structure of the hadrons. We start this investigation 
with the introduction of the entropy density s(T) for a thermal quark system at very low 
temperatures T. Here we assume the singlet structure for the confined quarks in the hadronic 
bag which should hold for temperatures well below the critical temperature T c . Herewith we 
are able to include some aspects of the groundstate structure in the equation of state. 

We start with the energy density s(T) and pressure p(T) of a confined colored quark gas 
for T <C T c . We use the First Law of Thermodynamics to include the classical heat density 
contribution 

s(T)T = e(T) + p(T) - fi q n q (T), (57) 

where n q is the quark chemical potential and n g (T) is quark density distribution function for 
a single quark flavor. We may rearrange this equation using the fact that the thermal average 
of the trace of the energy momentum tensor providess the equation of state as given in the 
equation (J2J), which now takes the form 

0£(T) = s(T)T - 4p(T) + N n q {T). (58) 

For the bag model in the limiting case that the temperature approaches to zero we use p = —B, 
where, as before, B is the bag constant [16] independent of the temperature. Clearly in this 
limit s(T)T vanishes. Furthermore, in the above equation of state we replace n q (T) with the 
Fermi distribution function f(fi q , T) times no, which is the quark number density in the bag at 
zero temperature. 

0£ o = AB + fi q nof(ji q ), (59) 

In the low temperature limit T) is simply a step function so that f(fJL q ) becomes the 

normalized particle distribution at zero temperature. As we have previously explained in 
the introduction, the trace of the energy momentum tensor can be related to the gluon and 
quark vacuum expectation values arising from the operator product expansion [16, 28, 36] 
following from equation © for 8£ . Thus we have, as before, included the important vacuum 
contributions of both operator dimensions three and four to the equation of state, which are 
both clearly independent of the temperature. 

In relation to the previous section on scalar meson dominance we choose as a simple 
special case that of the meson with the same light quarks surrounded by nuclear matter 

12 The author thanks Gerry Brown for his suggestions on the subject of the mesonic bound states and the 
problem of the meson mass [29] . 
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as an example for this investigation. Then we have from the Fermi statistics using the 
quark- ant iquark symmetry fiq = — (J> q and for the antiquark quantum density distribution 
function nq(T) = no(l — f(fi q ,T)) . Thus in the limit of zero temperature we have simply 
the expected value for the bag model [14, 52] 

% = 4B. (60) 

For the moment we take as a first approximation the case that we consider the value [16] of 
only the gluon condensate as about 1.95[GeV/ fm 3 ] . Here we use the fact that the values of 
the quark condensate for the light quarks are almost negligible. Thus we note that for this 
special case the bag constant B is rather big- just under 0.5GeV. However, if we take T 
finite and small, the gluon condensate does not change very much [28,37] until we reach rather 
high temperatures. 

Nevertheless, the constituent entropy at very low temperatures due to the color de- 
grees of freedom does have a very important contribution to the equation of state. In the 
singlet state it has been shown from the structure of the hadronic density matrix [97] that the 
quantum ground state entropy for the colored constituent quarks and antiquarks contribute 
an internal entropy with the value In 3. Although the hadrons themselves are pure states 
with zero quantum entropy, the constituents do have a finite entropy. In fact, in the singlet 
groundstate the quarks(antiquarks) have the maximum entropy of pure (color) mixing. Then 
the equation of state for the bag containing the constituents can be written with the additional 
contribution from the mixing of the colors in the quarks and antiquarks treated as separate 
consitiuents. The total constituent entropy density is then 2sm for both the quarks and the 
antiquarks treated as separate entities. 

Now we consider the total equation of state T£{T) containing both the groundstate 
and the thermal states using the decomposition in Equation ((TJ) together with the above 
Equations (JSHj) and (JBTJj) for the meson made up of the constituent quarks and antiquarks. The 
equation of state for the color singlet quarks in the mesonic bag is given by 

T;(T) = (2s M + s(T))T - 4(p(T) — B) + fj, q n (2f(n q , T) - 1). (61) 

At very low temperatures we would expect that s(T) and p(T) both to remain insignificant 
so that the main change in the equation of state is the value of 2sm = 1.824[l//m 3 ]. 
Physically we could relate this effect either to a lowered bag constant or a raised chemical 
potential. Thereby we may be able to see how much the finite temperature changes the actual 
value of the entropy due to the color mixing with its corresponding effect upon the equation 
of state. 

Along this same line we have also investigated a baryonic model at finite T includ- 
ing the ground state entropy for colored quarks. Here we have studied the relation of the 
temperature and chemical potential to B. In this previous work we calculated explicitly 
the effects of colored quark entropy on the bag pressure [98] in terms of the temperature. 
The general equation of state is rather similar to Equation (J61|) with the exception that the 
groundstate contains the quark entropy with the factor of 3s b- Following this work we looked 
into the effects of the quark chemical potential [99] on the entropy for a color superconductor. 
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VI. Conclusions, Deductions and Evaluations 



In this last part of our discussions of lattice calculations for the equation of state we 
bring together some of the main points in this report. An important aspect investigated in 
the Introduction was the breaking of both the scale and conformal symmetries when the 
strong nuclear interactions are present. This situation was very well stated by Stefan Pokorski 
in his book on Gauge Field Theories: "Scale invariance cannot be an exact symmetry of 
the real world. If it were, all the particles would have to be massless or their mass spectra 
continuous." 13 This statement in itself is not so very surprising when we are only considering 
massive particles like quarks or hadrons. However, because of the trace or conformal anomaly 
we saw in the second part from the simulation of pure gauge fields that this type of symmetry 
breaking continues for strongly interacting gluons even at very high temperatures. In the 
third part when the massive quarks are present in the simulations, the breaking of scale and 
conformal symmetries have a somewhat different relation to the phase transition because of 
the process of the restoration of the broken chiral symmetry which changes the properties of 
the quarks in the transition from the confined phase to the deconfined one. It is known, in 
fact, that the Wilson or Polyakov loop as the order parameter of the pure gauge deconfinement 
transition ceased to provide this property in the presence of dynamical quarks [2,3,5]. In 
this context we looked explicitly at the chiral anomaly in relation to the breaking of chiral 
symmetry in the presence of dynamical quarks. The thermodynamics is best related through 
the susceptibilies to the chiral condensate, which becomes itself an order parameter in the limit 
of vanishing quark masses. These results led the way to the explicit numerical evaluations of 
the thermodynamical properties of the gluon and quark condensates in the following sections 
of the next part. Finally we presented models which more directly relate to physical results 
from experiments. Now we summarize and expand upon some of the work which we have 
presented. 

VI. 1 Summary of Results and Related Ideas 

Throughout this work we have mentioned the importance of the ranges of the tempera- 
tures computed by using lattice gauge theory with the SU(N) color symmetry both without 
and with different numbers and masses of quarks present. By using the lattice data from 
numerous simulations we have plotted the various physical quantities in terms of the temper- 
ature. Without entering into the computational details we have investigated the numerical 
results in relation to the theoretical content of the simulation. We were in the above mentioned 
cases able to determine the thermodynamical functions like the pressure, energy density, 
entropy density and the resulting equation of state for the quarks and gluons as a function 
of the temperature. These results allowed us to see the growth of these thermodynamical 
quantities above the transition temperature, from which we could examine the quark and 
gluon condensates as well as compute the thermal properties in the determined related phases. 

There were a number of other important lattice studies of various different quark 
and gluon properties whose data and resuts we have not directly used. The starting point 
was in 1974 the "Confinement of quarks" by Kenneth Wilson [1,100] where he not only set 

13 This statement starts a "General Discussion" on the topic of Broken Scale Invariance [57] on page 237. A 
similar deduction was made on page 174 by Roman Jackiw [33] where he declared that "scale and conformal 
symmetries are broken, as they must be in order to avoid a physically absurd mass spectrum." 
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out the formalism of the lattice gauge theory by quantizing in Euclidean space-time, but also 
proposed the method of evaluation using the Feynman path integrals with the strong-coupling 
approximation. This very original work also introduced the Wilson action including the 
quarks 14 . An alternative approach [102] to the lattice fermions was proposed by John Kogut 
and Leonard Susskind known as the "staggered fermions." Most of the data presented above 
uses this formulation for the dynamical fermions on the lattice since it is somewhat better for 
the type of computations which have been carried out. The numerous works of Michael Creutz 
were very essential to these numerical simulations [2,103]. In the earlier times around 1980 
there were many very active groups which contributed greatly to the development of the lattice 
gauge computations with some very important results [104-107] for pure SU(N) lattice gauge 
simulations. These works provided the starting point for many of the later computations of 
the MILC and Bielefeld groups- for further literature see [11,12,103]. 

In recent years there have been considerable progress in the simulations with finite 
chemical potentials on the lattice. The above discussions have not included these many 
extensive computations [109,110] with dynamical quarks who included the chemical potentials. 
Only in the above very special example of the single flavored meson at low temperatures have 
we directly included the effects the quark-antiquark chemical potential. This particular case 
has its sole importance for the analysis of the model of a meson for the equation of state at very 
low temperatures. Others with dynamical quarks who included the chemical potentials were 
some works on the phase diagram by Fodor and Katz [109] as well as the more recent work 
done in collaboration of Bielefeld and Swansea [110]. These results we have not included with 
our analysis. There are also further important contributions from the CP-PACS Collaboration 
with two heavier quark flavors using the Wilson action for which we were unable to set 
the temperature scale to find thermodynamical functions out of the computed ratios [108]. 
Although all these works are very interesting, the form of the data as well as the values are 
often not in the temperature range where a comparison with the quantities are easily carried out. 

VI. 2 Implications from the Analysis 

There has oftentimes been a problem with the interpretation of some of the results 
from the numerous lattice computations due to the many general statements arising out 
of the ideal gas models. In numerous papers and even in some very highly respected text 
books 15 there have appeared statements where the expected approximate equality of e and 
3P becomes very nearly realized within a few integral mutiples of the critical temperature. 
However, we have seen above that within the range of the present numerical data from lattice 
gauge simulations all of statements of this type are clearly not fulfilled. In the earlier parts 
of this report we have seen the results for the equations of state in the figures [TJ |21 IH and 
All of these evaluations clearly show directly that these equations of state in terms of 
the temperature remain finite well above the critical or deconfinement temperature. As we 
have already explicitly pointed out in the text, with the possible exception of our figure El all 
the other plots show clearly a monotonical growth of the equation of state as a function of 
the temperature. Perhaps one possible source of misunderstanding of the lattice data lies in 

14 For a description of the actual Wilson fermions see the book by Jan Smit [101]. 

15 For example, at the end of a paragraph describing some numerical simulations [4] on page 10 it was written: 
"Even at T = 2T C , the energy density only reaches about 70% of the ideal gas value. However, at least for 
T > 2T C , one finds e w 3P." 
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the representation as ratios with the T N powers divided out. From the computational side 
this representation is very useful for direct comparisons with the high temperature limits for 
lattices of different sizes. However, for the interaction measure, which is defined as a ratio 
there are thermal contributions to the energy density which are different from those arising in 
the pressure. The difference e — 3p has contributions of lower powers and logarithmic terms 
in the temperature which appear to go to zero in the high temperature limit when divided 
by the higher power. Thus the interpretation coming from the appearance of the numerical 
results from the lattice presented as ratios can be quite misleading for the curves with a finite 
or vanishing limiting behavior. Here, as we have stated above, it is important that such terms 
be written in terms of the well understood thermodynamical quantites, which must diverge 
in the infinite temperature limit in accordance with the leading powers. Furthermore, as one 
perhaps could expect from the formulation of the laws of thermodynamics, we are also able to 
described the equation of state using exact differential forms which relate directly the known 
anomalies to the actual physical quantities [52,64]. In this general manner the action as a 
formal expression of quantum field theory is also clearly related to the general structure of the 
equation of state (see Appendix B). 

It was recently discussed how the lattice results show that even at temperatures 
much larger than the deconfinement temperature the thermodynamically related observables 
like the pressure ratio, the energy density ratio, the interaction measure and the baryon 
number density ratio at different chemical potentials still deviate by more than 20% from their 
respective relativistic ideal gas values [111]. These authors consider two special model cases: 
the two phase model and the mixed phase model, which they have quite appropriately named. 
The direct comparisons with the lattice data for the 2+1 flavor cases [70,71] and [109] provided 
insight into both the critical behavior and the asymptotic approach at high temperatures. In 
the case of the interaction measure the two phase model had a much sharper peak at the 
critical temperature than provided by the lattice data. However, the mixed phase model the 
peak provided much smaller peaks in the interaction measure. For the other ratios involving 
the pressure, energy density and the baryon number density the results of the comparisons 
with the lattice data were generally quite similar. 

VI. 3 Applications to Physical Processes 

We have already mentioned in the previous part of this report the properties of the 
equation of state at low temperatures for the study of the hadronic structure. In this 
temperature region the thermodynamical functions included in the equation of state are much 
harder to compute using the numerical simulations of finite temperature QCD on the lattice. 
Nevertheless, there has recently been considerable progress on the numerical computations 
of correlations between quarks and antiquarks at these lower temperatures in certain model 
cases [112]. Certainly a part of this computational difficulty lies in complicated quantum 
structure in the low temperature thermodynamics providing for the basic properties of the 
hadronic ground state. In this context we have previously looked into the role of the quantum 
entropy [97] in the discussion of the stability of the various quantum states- particularly 
relating to the singlet hadronic states. 

The thermodynamics of quarks and gluons in the confined region is of great interest 
in the examination of the formation of particle states below the deconfinement temperature. 
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The recent computations of the quark- ant iquark free energies [112] using the renormalized 
Polyakov loop has provided considerable new insight into the short distance interactions 
between the quarks and the antiquarks over a large range of temperatures. The results from 
these numerical simulations evaluate the color averaged and the color singlet free energies 
for static quark-antiquark sources placed in a thermal gluonic heat bath. An important 
procedure [112] in these calculations is the renormalization of these free energies using the 
short distance properties of the zero temperature heavy quark potential. This procedure leads 
to the definition of the renormalized Polyakov loop as an order parameter for the deconfinement 
phase transition of the SU(3) gauge invariant field theory. The color averaged free energy 
Fgg(r,T) of a static quark-antiquark pair acting as sources in a thermal medium is related 
up to a function of the temperature to the thermally averaged product of the color averaged 
Polyakov loops at two different spatial points separated by a distance r. The color singlet free 
energy Fi(r,T) is related to both the thermal and color averaged product of the gauge field 
variables. The significance of these results lies in the fact that at very short distances these 
two differently defined free energies become the same up to a factor of T In 9. Thus in the limit 
of short distances it was found [112] that 

hm[F 9q (r,T) - F 1 (r,T)} = Tin 9, VT. (62) 

r — >0 

However, at larger values of the distance r there are additional contributions to the colored 
averaged free energy coming from the octet states, which vanish at small distances. More 
recently the actual effective running coupling at finite temperatures has been calculated for 
both the quenched [115] and two flavor [116] QCD from the derivative of the free energy with 
respect to the separation r between the quarks and antiquarks. In a later work these latter 
authors [117] look at the singlet internal energy and entropy at large separations and higher 
temperatures about 1.3T C . However, in the case of small separations these authors note that 
both Lq(r, T) and TS\{r,T) suffer from lattice artifacts which result in the same problems for 
the singlet free energy. 

Next we mention some ideas closely related to these results but arise from quite dif- 
ferent considerations for low temperature quantum critical systems. First we recall from 
thermodynamics that the free energy F and the internal energy U are related through the 
entropy S by a simple transformation: F = U — TS, whereby we can write that S — — 
Thus we see that for finite temperatures we have a simple relation to the entropy. However, as 
mentioned in the last part, the entropy of the hadronic groundstate has the role of the mixing 
of the colors. It has been recently pointed out [113] that for one-dimensional quantum critical 
systems the entropy near criticality can be clearly divided into two separate parts- the bulk 
part arising from the partition function and the boundry terms relating to the groundstate 
degeneracy d g . This approach had been successfully applied in the past [114] to such important 
problems as the Kondo effect and the Heisenberg ferromagnet. Here we are interested in 
some different applications in relation to quark string models for the discussed computations 
on the lattice of Polyakov loop correlations [112]. We recall in our above discussion of the 
quantum ground state of the hadronic singlet state that the associated entropy of the quark 
and antiquark is In 3. If we regard the quark and antiquark as the end of a one-dimensional 
string, then we may associate the boundary entropy lnd g as arising from the groundstate 
degeneracy d g — 3 for the static quark and antiquark. Thus by the addition of entropies we 
may write the total boundry entropy for the two static sources as 2 In 3. All other sources of 
entropy we may associate with the bulk entropy in the string [113], which we assume here is 
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mainly due to the gluons. Finally if we assume that the internal energy at zero temperature 
is just the static quark-antiquark potential, we are thereby able to identify the term Tin 9 as 
arising from the groundstate entropy. 

We now recall some of the general results of our report as written above. The light 
quarks first go through the chiral restoration transition at a temperature T c then the effects 
of deconfinement actually will appear at a higher temperature The nature including even 
the order of these transitions appears to depend greatly upon the number of flavors as well as 
the masses. The heavier quarks show the separation between the transitions much less clearly. 
However, in all cases both the scaling and the conformal symmetries are still broken at much 
higher temperatures well above both of these transitions. We see from the physical equation of 
state that the very high temperatures provide a further symmetry breaking which is not present 
in the QCD vacuum contributions. The properties of the Brown-Rho scaling are important to 
the actual particle structure when surrounded by hot dense matter. Also recently the mass and 
the width of the sigma resonance have been calculated [118] by locating the pole in the pion- 
pion scattering amplitude with the quantum numbers of the vacuum 16 . The quarks and gluons 
are strongly correlated even well above the transition temperatures. There is also an effect of 
the running coupling [115,116] involving the confining part of the quark-antiquark interaction 
well above T c . In a very recent article Gerry Brown and collaborators posed a critical question 
on the results of the recent heavy ion experiments at Brookhaven National Laboratory. 
Therein they indicate the nature of the high temperature state by the statement [120]: "We 
suggest that the new form of matter found just above T c by the Relativistic Heavy Ion Collider 
is made up of tightly bound quark-antiquark pairs,..." Their observation gives further motiva- 
tion for the present investigations of the correlations and interactions at temperatures above T c . 

Finally in the two appendices we add some detail to the discussions of the models 
and the physical currents. There we formally discuss the actual physical meaning of these 
currents as expansion parameters and forces arising from the breaking of the scale and 
conformal symmetry. We take the well known example described in Appendix A - the MIT 
Bag model. The bag constant itself arises in the trace to make it nonzero. In the Appendix B 
we discuss how the dilatation and conformal currents are related to both the equation of state 
and the gluon condensates. We show these relations are best understood from the integral or 
the dual forms of these currents. 



16 This work has been further clarified by another even more recent investigation [119] of the sigma coupling 
to the photons in the process 77 — > tttt. 
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APPENDIX A: Phenomenological Models for Confinement 

In this added section we sketch the properties of some of the most prominent phe- 
nomenological models used in the theory of strong interactions to explain the observation of 
quark comfinement. It is meant as an extension to the introduction for the added properties. 

I. The MIT bag model was proposed as a model for the extended hadrons which con- 
tain free quarks inside a small volume. The action [13] with the internal energy density B is 
written in the following form: 

W = j^dt j fr[\tf - 1(V0) 2 - B], (63) 

where (pit, r) is the prototype of a hadronic field, which are the partonic or hadronic constituents, 
and R and (ti, t 2 ) are the space and time regions for the "bag". This action provides the field 
equations inside the bag 

-U<j>{x) = (64) 

and on the surface 

n-— 6 + n-V0= 0. (65) 
at 

Then the bag condition at the surface becomes 

\<? - \w? = B. (66) 

These three equations are basic to the MIT bag model. We note that the structure of this 
model is not directly in correspondence with QCD because of the time of development of the 
fields 17 . 

17 We shall use only the static properties of B for the special examples discussed in Appendix B 
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2. There have been many further developments on the equation of state for many cases 
involving bag models in a more general form. One can propose a form for the energy E(T, V) 
containing the bag term as E (V) = BV. As a correction one often adds the thermal radiation 
terms of the Stefan- Boltzmann type a SB T^V . However, this form may be extended to a more 
general energy equation of the following type: 

E(T,V) = E (V) + E'(T,V). (67) 

The pressure is defined by the usual relation p = —{dE/dV)T- When the free energy density 
f'(T) of the finite temperature part is scaled with a temperature d so that f'(T) = T<f)(i9/T), 
where is only a function of "&/T, then the equation of state is just 

p = -B + 7e'. (68) 

The constant is defined 7 = — d(ln i?)/d(ln V) for the energy density at finite temper- 
ature e'. For the ultrarelativistic gas 7 = 1/3, which yields the known relationship to 
the bag constant [13, 14]. The form of this equation for the relativistic gas is similar 
to the Debye equation of state for solids at low temperatures. In the Debye theory of 
solids 7 is often referred to as the "Griineisen constant." This approach of separating the 
groundstate and thermal states allows us to write a very general form of the equation of state 18 . 



3. The chiral bag model [21] has additional properties in relation to QCD. This new 
structure is often referred to as "The Cheshire Cat Mechanism." For simplicity one usually 
looks at a model in 1+1 dimensions. The action S is then decomposed into three parts 

S = S v + Sy + S dv (69) 

where for the fermions 

S v = [ d 2 x^d^ + ..., (70) 
v 



while for the bosons 



Sy = Jfo^d^) 2 + .... (71) 

The additional surface term becomes with the value of f — -k= and is the surface area 

J V47T 



with the normal vector so that 



~>dV 



[ dX'dnrfexptint/fW} + .... (72) 

JdV z 



APPENDIX B: Mathematical Forms and Physical Currents 

In this appendix we provide some additional information on some particular differential 
forms relating to the various discussed broken symmetries due to the properties of the strong 
interacions. There are two different types forms relating to the structures of the anomalies, 

18 The application of this equation of state to nuclear matter is described in the book of John Dirk 
Walecka [121]. 
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which are related to the axial and conformal anomalies respectively. In the former case there 
are already a large number of works [20, 34] relating to the chiral symmetry breaking and the 
axial anomaly, for which there is a very well known differential four-form AdQ. Furthermore, 
there is the related Chern-Simons Form dB, which can be represented as a three-form on a 
three dimensional closed surface. One can also discuss the physical interpretation of this form 
which is related to the topological charge. These ideas together with the properties of the 
chiral symmetry breaking in relation to the instanton solutions is discussed in the book of 
Edward Shuryak [15]. 

The scale and conformal symmetry breaking can be similarly related to the confor- 
mal or trace anomaly. In the above work we have found the presence of a finite trace of the 
energy-momentum tensor relating to the equation of state and the QCD sum rules. It is known 
in quantum field theory that classically the trace of the energy-momentum tensor is zero. In 
the following we derive (introduce) the corresponding three-forms under the names dyxle for 
the scale breaking and the fourspan for the special conformal symmetry breaking. The actual 
physical meanings for these currents and charges have been previously discussed [35]. The 
actual interpretion of these as energy flux and shearing forces has been given not so long ago 
by the author [64], upon which we shall now briefly expand. 

The dilatation current has already been defined above in terms of the position 
four- vector x M and the energy momentum tensor T^ u as simply the product xJT m as moments 
of the energy density. In the case of the general energy-momentum conservation [33, 35] one 
can find quite simply a relation to the equation of state [52,64]. We now look into a volume 
in four dimensional space-time Q containing all the quarks and gluons at a fixed temperature 
T in equilibrium. The flow equation holds when the energy momentum and all the (color) 
currents are strictly conserved over the closed bounding surface dQ of the properly oriented 
four- volume Q, which yields 



We have already introduced [64] the dyxle three-form as T>^dS^ on the closed three dimensional 
surface dQ, which is simply just the divergence in the four dimensional Minkowski space. The 
dyxle is the dual to the dilatation current written as the one-form D^dx^ in the usual four 
dimensional space-time 19 . In this context the dilatation current D^dx^ acts as if it were a 
streaching force over the space-time infinitesimal dx^. Whereas its dual form in space-time 
represents the divergence of this three dimensional closed surface acting as the boundry for the 
four dimensional volume fl On the right hand side of (|73*jl the integrated four-form j n T^dVL 
is an action or energy moment integral involving the equation of state. Since we assume 
here the positivity of the trace so that > 0, we can conclude that this action integral is a 
positive quantity. It acts as the color averaged source of the energy flux. This action integral 
gets quantized with the fields through the renormalization process. 

The analogous three-forms can be defined for the four special conformal currents 
which we shall collectively call the fourspan. The analogous dual forms are derived from the 
equation in a similar manner to that previously done for the dyxle [52,64]. The results for 

19 The relation of a differential form to its dual form is important here [122], see especially sections 4.6, 5.9 
and 10.6. 




(73) 
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the four conformal currents yield the four three-forms IC^dS^ is derived similarly from the 
following equation: 

I IC^dS" = [ 2x a T»d{l, (74) 
Jan Jn 

We point out here that for the fourspan the source terms are the first moments in space-time 
of the equation of state. Furthermore, the physical nature of the quantities /C° as the special 
conformal currents can be seen to be related to the shearing forces which destroy the conformal 
symmetry in all the four directions. The effect of the fourspan is analogous to the dyxle 
along the radial line in the four dimensions. In summary these five three-forms provide the 
mathematical structure for the known breaking of the scale and conformal symmetries for the 
strong (nuclear) interaction. 

As an example of the formal statements above we work out the analytical properties 
of these forms may be exactly calculated in the special case of the well known [13] MIT bag 
model. The scaling properties associated with the trace of its energy-momentum tensor we 
have noted above 20 . This relationship yields in the above equation (J2J) simply 9£ = 423. This 
result gives an exact solution for the dilatation current D^(x) of the form: 

D»(x) = Bx». (75) 

It is clear that this solution is a special case of the above general form x v T^ u for the bag model. 
Similarly we derive the special forms for the four conformal currents K afl (x) using the same 
general structure as for the dilatation current. Here there are two different types of solutions: 

(a) a = fi and (b) a ^ \i. In the first case (a) using the propertime r 2 = x^ with the implicit 
summation over the indices /i, we have for a single value of the index chosen from a — n — 
for the temporal case 

K 00 (x) = B(2(x ) 2 -r 2 ), (76) 
while for i — 1, 2, 3, where i = a = /i^0 the spatial case yields 

K u (x) = B(2(x 1 ) 2 + r 2 ). (77) 

For the case (b) we have t ^ k so that the result for l, k = 0, 1, 2, 3 is just 

K LK (x) = 2Bx L x K . (78) 

Thus we see that in case (a) in the equations ()76|) and we ()77|) have various combinations 
of quadratic terms depending upon which value of i — or i — 1, 2, 3. K 00 (x) is just a four 
dimensional shearing force. Furthermore, we notice that it is truly positive definite so that 
for both spacelike and timelike events in the bag it remains always positive. However, in case 

(b) in equation (fTSj) we find simply bilinear terms in two different space-time coordinates. 
The simplicity of these solutions for D^(x) and K a ^{x) arises form the fact that the energy 
momentum tensor and the metric tensor are both diagonal in the Minkowski metric for the 
bag model. All the conformal currents K atl (x) act like internal shearing forces which all go 
against the preservation of the angular symmetries which includes the breaking of the scale 
symmetry associated with D^(x). 



20 We have discussed in the Introduction the relation [14] of the trace to B as seen in their equation (1.10) 
therein and independently somewhat later [52]. 
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